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Abstract 



We study loop corrections to correlation functions of inflationary perturbations. Previous calculations 
have found that the two-point function can have a logarithmic running of the form log(A;//7,), where k is 
the wavenumber of the perturbation, and /j, is the renormalization scale. We highlight that this result 
would have profound consequences for both eternal inflation and the predictivity of standard inflation. 
We find a different result. We consider two sets of theories: one where the inflaton has a large cubic self- 
interaction and one where the inflaton interacts gravitationally with N massless spectator scalar fields. 
We find that there is a logarithmic running but of the form log(H/ fi), where H is the Hubble constant 
during inflation. We find this result in three independent ways: by performing the calculation with a 
sharp cutoff in frequency-momentum space, in dimensional regularization and by the simple procedure 
of making the loop integral dimensionless. For the simplest of our theories we explicitly renormalize the 
correlation function proving that the divergencies can be reabsorbed and that the correlation function 
for super-horizon modes does not depend on time (once the tadpole terms have been properly taken into 
account) . We prove the time-independence of the super-horizon correlation function in several additional 
ways: by doing the calculation of the correlation function at finite time using both the regularizations and 
by developing a formalism which expresses loop corrections directly in terms of renormalized quantities 
at each time. We find this last formalism particularly helpful to develop intuition which we then use 
to generalize our results to higher loops and different interactions. In particular we argue correlation 
functions have no long-term time dependence even if the spectator fields have a potential. 



The purpose of this paper is to compute loop corrections to inflationary observables. Why? At flrst there 
seems to be no good reason. Let us consider a standard inflaton (p with action: 



1 Introduction 




(1) 



In order to have an inflationary solution, all derivatives of the potential need to be small: 
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This means that the inflaton is extremely weakly coupled. In fact Maldacena showed that the leading 
interactions come from the mixing with gravity [1]. The size of the inflaton two-point function is of order 



{Sct>l)tree ~ ^ • (3) 

Oservations of the CMB imply that the amplitude of the primordial curvature perturbation is: 

where e is the slow-roll parameter and thus H/Mpi <C 10~^ during inflation. Loop corrections will be 
mediated by gravity, and will therefore be suppressed by Mpj. By dimensional analysis we expect: 

Wfe)l-loop~ -^-^ • (5) 

This is at least a factor of 10~^° smaller than the tree level result. Clearly nobody should bother 
computing it. 

This has been the general attitude for about the first twenty-five years after the invention of inflation 
until Weinberg stressed a different point of view. In the paper where he first studied these effects [2], 
he argued that since experiments seem to be providing more and more evidence in favor of the theory 
it is worth exploring all its predictions, even those that naively appear observationally unverifiable. 
Perturbation theory ought to be well defined even for inflationary fluctuations, and therefore, given the 
well known subtleties involved in studying de Sitter space, it is interesting to see how this works in 
practice. Many times as a result of doing a non-trivial calculation one can gain new insight into the 
theory and explore possible subtleties or generalization. For example, this was the case for Maldacena's 
calculation of inflationary 3-point function in standard slow-roll inflation. Though for standard slow roll 
inflation the effect was is in practice unobservable it opened the way to formulate alternative theories 
which predict much higher levels of non-gaussianities (see for example [3, 4, 5, 6, 7, 8, 9, 10, 11, 12]) and 
are already constrained by existing data [13, 14, 15]. 

The motivation articulated by Weinberg is of course a good one but we flnd two other issues to be 
even more important. The first is related to eternal infiation. This is a particular regime of infiation 
that occurs when the potential is flat enough so that quantum fluctuations dominate over the classical 
motion and there is a finite probability for inflation never to end. This regime of inflation has become 
more important recently because of Weinberg's proposal to explain anthropically the smallncss of the 
cosmological constant [17]. Eternal inflation provides a natural and simple mechanism, though not strictly 
necessary, to populate a number of vacua large enough to explain the tuning of the cosmological constant. 
The existent of such a phase was discussed in a more rigorous way only very recently [18, 19], following 
the original works in [16]. It was shown that if we start with a volume of space with the inflaton of (1) up 
in the potential there is sharp phase transition as we make the potential flatter and flatter: at a critical 
value of the slope, the probability of creating an infinite volume goes from being exactly equal to zero to 
non-zero. It was also shown that the number of e-foldings in any finite realization is bounded by S'ds/S, 
where Sds is the dc Sitter entropy at the end of infiation. This provides additional insight on how to 
interpret de Sitter space within quantum gravity [20, 19]. A fundamental ingredient of these proofs is 
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based on the fact that, at zeroth order in the slow roll parameters, the two-point function of the inflaton 
field at coincidence acquires a calculable logarithmic divergency which leads to a linear time dependence: 

/Aa[t] tt2 
d^k-^ ~ log(a) ~ Hh + const. , (6) 

where here A is a physical cutoff, and a{t) ~ e^* is the scale factor during inflation. Notice that the 
linear time dependence is in strict correspondence with the scale invariance of the two-point function. 
This linear time dependence was essential for the fact that eternal inflation could transition from being 
eternal to not eternal depending on the slope of the inflationary potential. If the time dependence were 
to be different, even by the slightest amount, inflation would be cither always eternal or never eternal. 
Loop corrections have the potential of changing the time dependence and is thus interesting to verify 
what the dependence is. 

Going beyond eternal inflation, it is well known that de Sitter space is a rather puzzling spacetime in 
the context of quantum gravity [20], and inflation offers a natural regularization for the inflnities of de 
Sitter space. Studying in detail inflationary spacetimes can lead to some new insight about its embedding 
into a theory of quantum gravity. 

The first calculation of the one- loop corrections was done in [2]. We find that result rather puzzling: 

(Cl)i-ioop = iV-^-^(c+gelogy) (7) 

where c is a numerical constant that depends on the renormalization procedure, e is the slow roll parameter 
e = —H/H'^, and /x is a renormalization scale. N represents the number of massless spectator scalar fields 
which interact with the inflaton gravitationally. It is this factor of N that parametrically distinguishes 
the effect from the spectator fields running in the loop from the effect of the gravitons. This offers a great 
simplification to the calculation ^ . The same form of the logarithmic running was found in all subsequent 
related studies [2, 21, 24] ^ 

Notice that if we take the above result at face value and we neglect the fact that it is slow roll 
suppressed then the two point function of the inflaton at coincidence would receive a correction of the 
form: 

/Aa{t) tt2 tt2 tt2 tt2 

d'A;^ • iV • ^ \og{k) ^-N-^.H^ {log{a)f ~ -iV • ^ • hH^ + const. 
PI ^ 

where we have neglected numerical constants but kept explicit the sign of the correction, which is 
negative. If we were to use this result to repeat the analysis of [18] we would find the very puzzling result 
that no model of slow roll inflation is eternal. 

Another puzzle with expressions (7) or (8) is that they seem to imply that for some large or small 
values of k or for some large enough time, the one-loop corrections become large enough to harm the 
perturbative expansion. This result seems to be at odds with the common intuition that, as a consequence 

^Here we are using the numerical result of [21] which corrected a mistake in [2] associated with the choice of 
the i e prescription. 

^It is worth to notice that usually in the just mentioned literature the above logarithmic running was simply 
stated as log(fc) -|- C, with C a constant containing the logarithm of a quantity with dimensions of mass. 
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of the approximate symmetry under time-translation and rescaling of the inflationary space time, every 
mode goes through the same history. The fact that the logarithm is proportional to a slow-roll parameter 
might suggest that it is precisely the deviation from de-Sitter which is generating the effect. A closer look 
at the calculation reveals that this is not the case. Furthermore such an effect would imply the existence 
of a non-trivial correlation between very different modes which we find hard to believe. 

An obvious problems with the result in (7) is that the argument of the logarithm is a ratio of a 
comoving scale k and a physical scale /i. This cannot be the case and it is the sign of a mistake in the 
calculation. In fact by its very definition the scale factor of an FRW Universe is unobservable. It just 
corresponds to a rescaling of the spatial coordinates. This implies that every physically correct result 
should be invariant under the following rescaling by a real number A: 



If we express the result of (7) in real space, we have to do two momentum integrations by d^k d^k' (in 
(7) we neglected the momentum (5-function which of course is there because of translation invariance). 
At this point the prefactor of the \og{k/ jx) is invariant under (9) while the \og{k/ jj) is not. 

The minimal way the above expression could be corrected is by replacing the argument of the log 
with k/{ait)ij). But what value of i should we use? There are basically three options: i ~ t^arly^ i ^ t 
OT t ^ tk, where tgariy is some early time before the mode crosses the horizon, t is the time at which the 
correlation function is being calculated and tk is the time when the mode k crossed the horizon. The first 
option is ruled out by adiabaticity: the mode in the past has a very high frequency and is in its adiabatic 
vacuum. It therefore has no memory of the early times. The second option would imply an additional 
time dependence in the two-point function which would result in no model of slow roll inflation being 
eternal. Furthermore after a very long time (but no time is too a long in the case of eternal inflation!) it 
would also mean that the one loop corrections would become large, comparable to the tree level result. 
This is at odds with the common intuition of an approximate time-translation and rescaling symmetry 
of the inflationary space-time as well as with the expectation that very different modes should be only 
trivially correlated. More importantly as we will soon see if this were to be the case it would create 
problems not only for eternal inflation but for inflation in general! 

If instead the true result is the third one then since k/a{tk) ~ H the logarithm would just become 
\og{H/ jj), a result that would make a lot of physical sense. This can be understood by recollecting how 
logarithms appear in scattering amplitudes. They appear in the form of log{E/iJ,) where now E is the 
invariant energy of the collision. This logarithm takes into account the difference between the theory 
defined at the renormalization scale /j, and at the scattering energy E. The logarithmic corrections from 
the one-loop calculation are not large once the renormalization scale is chosen to be close to the scale at 
which the experiment is done, which led to the invention of the renormalization group [25, 26]. We can 
reach a similar conclusion here: the energies probed by the interactions during inflation are of order H 
and therefore the logarithms are small if we renormalize the theory at that scale. 

Before moving on to the actual computation, we would like to stress a second motivation for studying 
loop corrections that we also flnd very important. There is another disturbing consequence that would 
arise if the way of making the result in (7) correct was to replace the argument of the logarithm with 
k/{a{t)ii). The theory studied by Weinberg and giving rise to (7) is a theory of the form 



a — > A a , 



X ^ x/X 



k 



Xk . 



(9) 




(10) 
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where there are massless scalar fields ai interacting gravitationally with the inflaton, sourcing the 
fluctuations of the metric C, through the vacuum fluctuations of their stress tensor T(^a)iJ.u- Since these 
spectator fields are massless, only their gradients are important and therefore the vacuum fluctuations 
of T(^a)iMi/ s-re expected to become uncorrelated for distances longer than the horizon scale. If the correct 
answer were log(fc/(a(t)/i)) it would mean that these fluctuations are capable of sourcing a fluctuation 
even on scales much longer than the horizon. This would be very surprising and would create serious 
problems for the theory of inflation. In fact the prcdictivity of inflation relies on the fact that the 
perturbation C is expected to be constant on scales much longer than the horizon independently of what 
happens on scales of order of the horizon. This is very important because for some epochs such as 
reheating or a GUT phase transition (if this exist), we have almost no idea of what happens on scales of 
the order or shorter than the horizon. In principle large fluctuations on horizon scales can exist during 
these epochs. They are correlated only over horizon scale distances and so one expects they cannot affect 
a mode of much longer wavelength. Instead a time dependence from the one loop corrections would 
imply the contrary. Notice that the situation can become really worrisome. The effect of the a fields, even 
assuming a time dependence, would be very small in standard inflationary scenarios. This is so because 
there is not enough time to overcome the tremendous suppression of order H'^ /Mp^. Basically the a flelds 
are free flelds and therefore do not have large fluctuations. The situation might be completely different for 
fluctuations at the time of reheating or at a GUT phase transition where the relative fluctuations might 
well be of order one even on Hubble scales. This could induce a change of ^ of order one in a Hubble 
time: a huge effect that would undermine our capability of making predictions out of the inflationary 
phase without having a detailed understanding of all the epochs in between us and reheating. The same 
effects would probably also alter the scale invariance of the primordial density perturbations. 

The purpose of this paper is to show that indeed there is no time-dependence nor deviation from scale- 
invariancc induced by these interactions and that the one-loop logarithmic correction to the two-point 
function is of the form \og{H/ ji). Achieving this will be a rather challenging and sometimes technical 
task and for this reason we will give several independent derivations, which will also help us in building 
intuition. 

After reviewing in sec. 2 the formalism for computing loop corrections, in sec. 3 we will start by 
studying a different and simpler theory given by [3] 



(11) 



where tt is the Goldstone boson of time translation, that are broken during inflation, and is related to 
by the simple relationship ( = —Hit. We will comment later on how this theory is derived. Here we 
just would like to point out that the scale M which controls the size of the self-interactions can make 
them parametrically larger than gravitational. In fact, its current upper limit is just observational: the 
cubic self-interaction induces a three-point function in the CMB and is thus constrained by WMAP [14]. 
Therefore we can study loop corrections induced by this self-interactions without including gravity. We 
think this is technically a simpler example. 

We will calculate loop corrections to the two-point function in this model. We will perform the 
calculation with two kinds of regularization. The flrst, in sec. 3.1, will be imposing a flxed physical 
cutoff in frequency and momentum. We will find some power law divergencies, finite terms, and most 
importantly logarithmic divergencies. As usual, the form of the logarithmic divergencies is computable 
in the infrared even without performing an explicit renormalization, and it represents the running of the 
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inflationary two-point function. After we remove the power law divergencies through renormalization 
the logarithmic corrections can become parametrically much larger than the remaining finite terms if the 
renormalization scale is taken to be very different than the energy typical of the inflationary processes 
which is Hubble, up to the point of potentially harming the perturbative expansion of weakly coupled 
theories. This is a situation historically known as the 'large logarithms'. Then, in sec. 3.2, we will perform 
the same calculation in dimensional regularization (dim. reg.) finding the same result for the logarithms. 

Dimensional regularization was the scheme used in all the literature that followed Weinberg's first 
paper. We point out explicitly where we believe a mistake was made and perform the correct calculation. 
Still working in this regularization, in sec. 3.2.1, we will be able to show by making the loop integral 
dimensionless that even without explicitly doing the calculation one concludes that if a logarithm is 
present it has to be of the form log{H/fi). Finally, in sec. 3.3, we will perform the renormalization of the 
two point function by identifying and solving for the proper local counterterms. 

In this process we will be able to address another issue that arises in these one-loop calculations. 
Because of the symmetries of the space time, we will perform the loop phase space integrals using a 
partial Fourier basis where we keep time real but use spatial momentum. In doing these integrals it has 
been noticed that because of the nature of the interactions, the integrals in time converge even when the 
external time is taken to infinity. It is only the integrals in momentum that have UV divergencies, which 
are supposed to be reabsorbed by local counterterms. The calculation in dim. reg. becomes simpler if the 
integrals are done in this order. However, as it was noticed by Weinberg [2, 22], the momentum integrals 
that are left can be more divergent than the ones one would get if the time integral was done up to finite 
time or if the order of integration had been switched. This implies that there is no guarantee that all 
the divergencies, including the logarithmic ones, obtained with the outlined procedure are reabsorbable 
by local counterterms. It is in fact possible that they would disappear if one where to switch the order 
of integration at the cost of inducing a time dependence of the one-loop result. This issue had been 
simply mentioned but not solved in the literature. We are able to address and solve this issue in several 
ways. Of course, the more direct way is by explicitly performing the renormalization of the correlation 
function, as we do in sec. 3.3. The second way is by performing the calculation regularizing the integrals 
with a cutoff, in which case we first integrate in momentum and then in time finding the same result as 
in dimensional regularization. This means that the divergency that we find has to be reabsorbed by a 
counterterm. Third in sec. 3.2.2 we actually perform the dim. reg. calculation keeping the external time 
finite finding again the same result. 

At this point extending our calculation to the theory of (10) is quite simple. We do it in sec. 4, 
finding again a logarithm of the form \og{H/ jjb). For this theory, in sec. 5 we prove that there cannot 
be any time dependence of the two point function in yet a different way that we also then generalize to 
the case of the vr self-interactions. This last proof can be thought of as computing the loop corrections 
by first doing the momentum integrals and then doing the time integral. In this way, we can assume 
that all the UV divergencies have been reabsorbed by a counterterm. This is so because in this way of 
doing the calculation, we are using directly the physical (renormalized) quantities defined at each time. 
In this method, computing loop corrections becomes equivalent to solving perturbatively the equations 
of motion for the various operators. It turns out this alternative method is not powerful enough for 
us to determine the exact result of the loop calculation but it allows us to show in a very simple and 
physical way that the C, two-point function is time-independent. In particular, the reader not interested 
in the technical calculations but on the main result could read directly sec. 5 to see that the C, correlation 
function is time- independent and then make the integrals dimensionless as in sec. 3.2.1 to see that the 
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logarithm, if it exists, is of the form \og{H/ ji). 

We find this last method of dealing with loop corrections particularly physical and intuitive. In 
section 6, we use this intuition to argue how our results extend to higher loops and different interactions. 
In particular we will claim that contrary to what was reported in the literature [22, 27] the Q two point 
function does not acquire any time dependence even in the case where the spectator a fields have any 
form of stable potential ^. It is conceivable that our results can be extended to other time dependent 
effects found in the literature when studying the zero mode in different theories [28] . We will summarize 
our results in sec. 7. 

2 One-loop corrections 

Let us write the metric in the standard ADM parametrization: 

ds^ = -N^dt^ + hij {dx' + N'dt) {dx^ + N^t) . (12) 

We can fix the gauge by choosing the inflaton (f) to be uniform on equal time slices: S(f) = and by 
imposing the spatial metric hij to be of the form: 

hij = a{tfe^^5ij. (13) 

In inflation we are interested in computing late time expectation values of the metric fluctuation C,. 
This is due to the fact that we expect that for wavenumbers much outside of the horizon it is constant 
in time and more importantly independent of the local, often unknown, physical processes that happens 
during the various epochs of the history of the Universe. Here we are going to compute the one-loop 
corrections to its two-point function. This includes checking for its constancy in time. The expectation 
value for C, is given by [1, 2, 21] 

{C\t)) = {Uintit, -00+)Wt)Uint{t, -00+)) , (14) 

where 

C/,„,(t,-oo+) = Te-'^-+'^*'^-*(*\ (15) 

where Hint is the interaction Hamiltonian in the interaction picture, T denotes time-ordering, and — 00+ 
means that the integral is performed on an analytically rotated contour t' — > t'{l + ie), which projects 
the free vacuum state onto the interacting vacuum state in the infinite past. The operator on the 
right-hand side is meant to be the freely evolving operator in the Heisemberg picture. At one-loop, 
eq. (14) becomes 

ft rt2 

(C'(O)i-ioop = -2Re' ' 



/ dt2[' dt,{H^^{t,)H^,{t2)e{t)) -2Im / dt,{H^,{t,)e{t)) 

r dh f dt^{H^,{t,)em^,{t2)) , (16) 

J—OO- J — 00+ 



■^This of course includes and agrees with the main conclusions of [22, 27] about the absence of an exponential 
dependence in time in the two-point function. 
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where -f^f^] and H^^^ represent the cubic and the quartic interaction Hamiltonian. We notice that there 
are two terms that involve two insertion of a cubic interaction, and one term that involves one insertion 
of a quartic interaction. 

We are going to apply eq. (16) to two different physical Lagrangians that represent two different 
models of inflation. The first is a model of single field infiation where the inflaton is self-interacting with 
a strength that is parametrically larger than gravitational, and the second is a model where the inflaton 
is gravitationally interacting with N free massless scalar fields. The fact that there arc N identical scalar 
fields makes the effect coming from them as they run in the loop parametrically larger than the analogous 
effect from loops of gravitons. 

Before actually proceeding, let us specify what is our target for the calculation. The one-loop correc- 
tions in eq. (16) will contain divergent pieces and finite pieces. It is expected that the divergent pieces 
will be reabsorbed by counterterms in the Lagrangian, which will affect also the finite terms of the result. 
In the present paper we will not perform the explicit renormalization of the two Lagrangians we will 
study, but we will rather assume that this can be done ^. For one kind of divergencies, the logarithmic 
ones, it is particularly easy to realize how they will be transformed by the renormalization procedure. 
This is so because a logarithmic divergency will appear in the form 

^'°K somephysiotl.na.s scale ) ' <^^) 

where C is the coefficient of the logarithmic divergency, A is the cutoff of the calculation, and in the 
denominator we have some physical mass scale. Upon renormalization, the above expression will simply 
be transformed by replacing the cutoff A with the renormalization scale /j,: 



A 



^■°H 3omaphysMn,„.sscde j ^ ^ '"H soma physM n,„.s scale j ' 

This means that even without explicitly renormalizing the theory, we can compute the coefficient C of 
the logarithmic running. Notice however that in principle we are unable to compute precisely what the 
mass scale which accompanies the renormalization scale is as we can always redefine it by a multiplicative 
constant at the cost of simply changing the finite terms. These finite terms are fixed once and for all 
by imposing the two-point function to have a certain value at some energy scale. They arc part of the 
definition of the theory. The logarithmic divergency tells us instead how the amplitude of the process 
changes as we change the renormalization scale: it represents the running, in this case, of the two-point 
function and it can be large if the renormalization scale is chosen to be very different from the relevant 
energy scale of the process. We will argue that there is only one physical mass scale that can possible 
accompany the renormalization scale in this situation, while at the same time leaving the remaining 
numerical finite terms of order one. Several calculations of the coefficient C and of the mass scale that 
accompanies the renormalization scale have been done since Weinberg in [2] , finding that the form of the 
log divergency is 

k' 



Clog(-) , (19) 



■^In particular, for one of the theories, we will explicitly verify that this is the case by identifying and computing 
the counterterms. 
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where k is the comoving momentum of the modes of which we compute the expectation value. We find 
this result not to be physical, as the comoving wavenumber k is not a physical quantity. Its value can 
be changed by simply a constant rescaling of the scale factor a, which is not physical. As we pointed 
out in the introduction, this result could have had disastrous consequences for slow-roll eternal inflation, 
and more in general for inflation as well. The main point of this paper will be to show that all of those 
calculation are not correct, and that the mass scale accompanying the renormalization scale is the Hubble 
scale H: 

<7 1og(^) ^ Clog(^), (20) 
proving in the same time that there cannot be any time dependence in the correlation function. 



3 Large tt^ self-interaction 



In [3] an effective field theory for infiation was developed that made it possible to explore in full generality 
all the possible self-interactions of the infiaton, and in particular to clearly see that the infiaton can have 
large self-interactions without spoiling the background quasi de-Sitter solution. It was shown that the 
self-interactions can be parametrically larger than the one mediated by gravity and therefore we can 
concentrate on those without worrying about metric fiuctuations. This simplifies the calculation a great 
deal. The Effective Lagrangian is a function of the field n which represents the Goldstone boson of time 
translations which are spontaneously broken during inflation. It is related to the ( metric fluctuation at 
linear level by the simple relation 

C = -Hn . (21) 

Notice that tt is a scalar held with dimensions of inverse mass. Here we do not present the derivation 
of this Lagrangian, that can be found in [3], and we simply take a limit of its parameter space and use 
it for computing loop corrections. There are in general two interactions of comparable strength, tt^ and 
7r(5i7r)^, each one accompanied by quartic interactions due to symmetry reasons. Upon tuning [3, 14], 
we can make the term containing vr^ parametrically larger then the one containing 7r(3i7r)^. This allows 
us to reduce the number of vertices, and therefore to simplify the calculation. In this paper we are trying 
to assess conceptual and qualitative features of the loop corrections to the inflaton two point function 
and we can therefore study a tuned theory. Following the notation of [23] , the action is of the form 



S= (fx 



-HM, 



PI 



1 



TT 



(5.7r) 



27r^ + 



(22) 



Here C3 is a dimensionless parameter, while M has dimensions of mass. The interaction picture Hamil- 
tonian is of the form 



int 



(f'x 



1 



PI/ 



(23) 



We notice that we have a very simple interaction Hamiltonian at cubic level: only one interaction of the 
form TT^. 

We are now ready to begin the computation. We will perform it in two different regularization 
schemes: the flrst is by putting a sharp cutoff in frequency and momentum space, and the second by 
dimensional regularization. 
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3.1 Regularization by Frequency and Momentum Cutoff 

(3) 

We start by computing the two contributions with two insertions of H-^l_, and in particular the one of 
the form: 

r /•*2 



{Cit))i-ioop,A = -2H'Re 



[ dtj' dt,{H^,{t,)H^,{t^)^\t)) 

J—OO- J —OO- 



that we refer to as contribution A. After some simple algebra we are left with: 



(24) 
(25) 



X 128 ciM^H^ Re 



t2 



dt2 a{t2r I dh a{hf I [ ^(2;r)3<5(3) (k + h + k) 

■cl/j. -cl/j. -cl ( J. \-cl*(j. \ -cl /j. \-cl*i 



X ^^'(ti)^r(i) ^t\t2)i:Tit) ^t[{h)nt[*{t2) <(ti)<*(t2: 

Notice that there are other possible contractions of the operators, giving rise to diagrams that are not 
one-particle irreducible. However they involve derivatives of the wavefunctions at A; = which vanish. In 
order to get the coefficient of the one-loop logarithmic divergency, we realize that we can simply regulate 
the above integral with a cutoff in frequency and momentum space. Though this regulator is not diff. 
invariant we can exploit the fact that the logarithmic divergency receives equal contributions from all 
energy scales and therefore its coefficient is independent of the UV regulator. The coefficient of the 
logarithm obtained in this way is the correct one at one-loop but not at higher loops as in this case it 
depends on the finite terms of the lower loop calculation which in turn depends on the regularization and 
renormalization procedure. 

We therefore regulate the momentum integral with the following replacement: 

d^k / d^k , (26) 

where A is a fixed physical cutoff, and therefore it appears cutting off the integral in comoving momentum 

k accompanied with a factor of the scale factor a, so that /sphysical = k/a{t) < A. Notice that the scale 

factor a is evaluated at time ti. This can be understood in the following way. If we look back at 

eq. (14), we see that it involves the interaction picture evolution operator U acting on the vacuum. In 

the interaction picture, we can think as the operators evolving with the free Hamiltonian, while the state 

evolving with the interaction Hamiltonian. The diagram we are computing corresponds to letting the 

vacuum state on the left (and then also the one on the right once we take the real part) evolve up to time 

(3) 

t with two insertions of the interaction Hamiltionian the first at time ti and the second at time t2, 
with ti <t2- Momenta which are higher than the cutoff at the earlier time ti, do not contribute to the 
vertex at t\ , even in the case that by the time t2 they have redshifted inside the cutoff. In the case there 
is no interaction at ti, then there can no be a subsequent interaction at time t2 and therefore we have to 
put the cutoff as Aa(ii). We have also to regulate the integral in ti. This is analogously regulated with 
the same physical cutoff by replacing the ti integral in the following way: 



/t2 



/ 



dU 



(27) 



At this point, the calculation becomes quite straightforward. Before actually doing it let us first gain 
some intuition by doing the calculation in a non-expanding and then in a slowly expanding Universe. 
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• Non-expanding Universe 

We start by completely neglecting the expansion of the Universe and take the wavefunctions to be: 

where we have set the scale factor a to be equal to one. Notice that H enters here just to canonically 
normalize the field, and does not play any role in our discussion. A straightforward calculation gives: 

(C,-(i)C,-,(0)l-loop, A,No Expansion ~ {27r)'S(^Hk + ^O^^^p^' (log (f ) + , (29) 

where L is a divergent constant and where we have neglected numerical factors. The result is clearly not 
scale invariant. Since the vr'^ interaction is of dimension six, the only way it gives rise to a logarithmic 
divergency is by pulling out some powers of the external momentum k. Notice that as A; — > 0, the one-loop 
correction goes to zero, as the interaction we are dealing with has at least one derivative acting on the 
external leg. 

• First corrections from the expanding Universe: kr] ^ 1 

Let us now include the first corrections from the space expansion and use the wavefunctions in de-Sitter 
space, taken in their high energy limit fc?? 3> 1, where 77 is conformal time: 

where in the second passage we have taken the krj ^ 1 limit. In this case the expansion is very slow for 
the modes being considered. The result of the integration gives: 

(CE{t)4,(*)>.-,o„P,A,.,». ~ (2n)V^Hk^k')^^^ (log +C 

(31) 



We see that, apart for the 1/k^ term and the (5-function, every other k is accompanied by an 77. In this 
way the expression, once Fourier transformed to real space, is correctly invariant under the rescaling 
a ^ Xa, k ^ Xk, x ^ x/X, which is a symmetry of the problem. Notice that again, in the limit kr] 0, 
the one-loop correction goes to zero. This is due to the fact that in order for this operator to be able to 
produce a logarithmic divergency, some powers of the external momentum have to be pulled out of the 
integral. Symmetry arguments dictate that each additional power of momentum must be accompanied 
by a factor of 77. This forces the expression, in the limit 77 ^ 0, to go to zero. Further the structure of 
the logarithmic divergency has changed: now the factor of k is accompanied by a scale factor evaluated 
at time r). We are beginning to see the relevant effect of the expanding Universe. We can anticipate 
what the effect of including the expansion of the Universe even at later times {krj ~ — 1) will be. As 
the expansion of the Universe goes on, and kr) becomes more and more infrared and eventually of order 
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one, the amplitudes of the modes freeze and are evaluated when kij ~ —1. This implies that the result 
becomes scale invariant, and that the logarithmic divergency takes the form of 

log(f) . (32) 

De-Sitter space has provided an infrared cutoff to the loop corrections. At the tree-level, quantum effects 
on scales larger than 1/H are larger than what they would be for the same scales in Minkowsky space. 
The same happens here for the loop corrections: on large scale, quantum effects are parametrically larger 
then what they would be in Minkowsky space. 

Notice two further comments. First, at least in principle, the term in (32) is a local term, and therefore 
one could imagine that a local counterterm proportional to log(i^/ iJ.) could remove it completely. However 
this is not the case because for modes that are still inside the horizon the divergency is of the form of 
log(k / (a(t) fi)) which is not a local term and therefore cannot be removed by a local counterterm. The 
counterterm should be the same for all the modes, and therefore we conclude the log(H/ jj) is something 
that cannot be removed by a local counterterm. Second, the way the logarithm in (31) becomes (32) as 
the mode goes outside of the horizon tells us that the H in the logarithm should be interpreted as the 
energy at which the process is evaluated for modes that have exited the horizon. This implies that we 
expect the constant C to be of order one once ji is chosen to be of order H. 



• Calculation in the Inflationciry space 

We are now ready to do the calculation in an inflationary Universe. We use the following classical 
wavefunction, expressed in conformal time: 

We would like to stress the technical point that since the momentum cutoff is time-dependent once 
expressed in comoving coordinates one has to perform the momentum integrals first and the time integrals 
second. As we will see in the next subsection, in dimensional regularization we will perform the integrals 
in the opposite order. Taking the t — -|-oo limit a straightforward calculation leads to: 

(Cg(t)C,-,(i))i-ioop,A, w+oo = -(2-)'^^'^(^>^^')p^^^log [j) , (34) 

where we have neglected all the finite terms and the power divergencies. Some details of this calculation 
are given in App. A. The result is time-independent and scale invariant. The other contributions do not 
produce any other logarithmic divergency and we will be left with log(iJ/A). 

Let us now proceed with the remaining terms. Let us analyze the additional one involving two 
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insertions of H^^^ which we wiU label with a B: 

(Cfc(i)Cfc,(i))i-ioop,B = i^' r dh f dt2{H^t{U)'K^{t)'K^,m'£S2)) (35) 

J —oo- J —00+ 

X 12S clM^H^ f dhaihff dt^ a{t^f j ^ J ^^{2^f 5^^^ (k + k + k^) 

X i^tit.W^t) 7rt,*{h)7rt,{t) <(tiX*(t2) <(tiX*(t2) . 
This integral is convergent and gives 

1 SSI r^H^M^ 

(C,-WC,-,W)i-ioop,B, t-.+oo = (27r)V3)(g^fe.^ 3 . (36) 



A;3 72007r2 H'^M^^ 



Finally, we have to compute the contribution from the quartic interactions. It is quite easy to see 
that they do not give rise to logarithmic divergencies. Here for brevity we just concentrate on the tt^ 
interaction, as this is the only one that will give an effect proportional to c| and so it could affect the 
logarithmic divergency we found in the first diagram. Keeping only the term proportional to C3, the 
expression, that we call contribution C, reads: 



/- t -| 2r7-2?i^8 



X Im 



ft pAaiti) j3h. 



PI 

(37) 



A straightforward integration shows that this contribution contains no logarithmic divergencies and goes 
as 

(C^(i)C,-,(t))i-ioop,c, t-.+oo = i2nf S('Hk + k')^sj^O (A^ H'A') . (38) 



Summarizing, we conclude that the two-point function receives a logarithtmic correction of the form 
log(i7/A): 

(C,-(i)C,-,(*))i-ioop, = {2nf6^'\k + ^^Op^^i^ (log [j) + Cj , (39) 

where C is a numerical constant. We note that the unitarity bound Au of this theory is of order 

which means that the above one-loop correction scales as 

(C')l-loop~(Otree(£)' - (41) 

As expected, the theory becomes strongly coupled when H ~ A[/. 
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3.2 Dimensional Regularization 

We now perform the same calculation as in the former subsection but in dimensional regularization. 
Dimensional regularization has the advantage of respecting diff. invariance. As we said before, this 
is not really necessary for computing the coefficient of the logarithmic running but it would be a great 
simplification if one wished to renormalize the theory and compute the finite terms of the loop corrections. 
Further, this is the regularization scheme that has been applied in the literature on the subject [22, 21, 24], 
and we claim that it has not been apphed correctly because several terms coming from the generalization 
of the expressions to d dimensions have been omitted ^. For these reasons we proceed to perform the 
same calculation in dimensional regularization (dim. reg.) where we make the integral convergent by 
changing the number of spatial dimensions. 

Let us start with the diagram A, which becomes 

(Cfe(i)C,- (i))i-ioop,A = -{27rf5^^+'\k + k') USclM^H^fi^' (42) 

X i^t{ti>k*{t) ^t'{t2)^k*{t) <(iiX*(i2) <(iiX*(t2) 

where 5 represents the difference between the number of spatial dimensions d and 3 and ji is the renor- 

malization scale and it has been inserted in order to keep C, dimensionlcss and vr with dimensions of time. 
Here we do not extend to d-dimensions the numerical factors as they would only change the constant 
terms from the loop-corrections. At this point we have to express the wavefunction tt'^' in d-dimensions. 
This reads: 



^fc / . \l/2 t.(3+5)/2 

r(i) 



ki3+S)/2 ^(3+5)/2( M ' (^^) 



where Hi> ' is the Hankel function of the first kind of index u. We do not actually need to perform 
the integral in (42) with this wavefunction. In fact, the two wavefunctions which just depend on the 
external time t and the external momentum k can be brought outside of the integrals and taken to three 
dimensions and to the ?7 — > limit. They give: 



i e^^^l^ H^+^l^ 1 i H 1 

^ 0) - ^^^^^^2 (_ijM2j/i<5)V2 A;(3+'5)/2 ^ (27r)V2 (-HM^^y/^ kV^ 



This is justified by the fact that every logarithm coming from the external wavefunctions will be canceled 
by the contribution from the counterterm which renormalizes this interaction. The same is true for the 
(5-function of the external momenta: 

S^^+^\k + k') ^ S^^\k + k') . (45) 

For the remaining wavefunctions inside the integrals we notice that we are interested in computing just the 
logarithmic running from the loop integral and we can therefore Taylor expand the above wavefunction 



^We will see that having forgotten those terms is equivalent to putting a cutoff in physical frequency-momentum 
space that grows exponentially in time. 
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around 6^0, and keep only the linear term (as we will see, the loop integral has only one simple pole). 
Notice that the momentum in the wavefunction is still the full d-dimensional one which is a sufficient 
condition for keeping the integrals regularized. Taylor expanding (43) around (5 — ^ we obtain: 

-t' - V . ,y2 ^2(^ + ^'^v)e-^'^ X (46) 

(l + . (i log i-H,) + ^ - l^e-^ (-3.. + 2 Ci(2.,) - 2 ^ 81(2^)) 

(-2iJM2j/x'5)V2 ' 

(^1 + 5 Q log {-Hr]) - ^e^'^'^ {-Zm + 2 Ci(2A;r?) - 2 i Si(2A;r;)) + O ((5^) 

Ci(x) = r dx' ~ ^ + log(x) + 7 , Si(a;) = f dx' 

Jq x' Jq x' 



where 



(47) 



and 7 is the Eulero-Mascheroni constant. We can now begin to do the integrals in (42) by starting from 
the ones in time. 

• Contributions at order 6^ 

We start by performing the time integrals in (42) for the component which is explicitly independent of 
6 (notice that k is still the d-dimensional momentum). This was the only contribution included in the 
literature. We get: 

q „2 zt8 i\/r8 

(C,-(t)Cg,(t))i-ioop,A,t^+oo,.o = -i^^fS^'\k + k')—^A—— (48) 



167r3 H'^Ml.^ 



1 f d^+^ki f d^+^h 

X 



^ ^ X kik2 U {kl + k2f + 9{ki + k2) k + SkA 

(2.)3«*<3«) + fa + fa) ^ 1- , 



fe7 J (27r)3^5 J (27r)3 ^ ' V ' V ^k + ki + k2 

where the subscript 5^ reminds us that we are taking only the terms in the wavefunctions of order 5^ . 
Notice that in this case, contrary to what wc did in the case of a regularization with a sharp cutoff in 
frequency and momentum space, wc perform the time integral up to t ^ +oo first, and then we perform 
the momentum integral. We will come back later in sec. 3.2.2 and sec. 5 to prove that there are no 
subtleties associated with this. By dimensional analysis, the remaining momentum integral will give a 
result of the form 

fc^f^VF(<5), (49) 



Because of the ultraviolet divergencies F{8) has a singularities in the limit (5 — of the form: 

F{5) = ^ + F^. (50) 
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In the limit 5 ^ it leads to: 

{logik/ fi) + L) , (51) 
where L is a divergent constant. It is straightforward to carry out the momentum integral and find: 

(C,K*)C,-,W)i-ioop.A,w+oo,.o = -{2.f5^){k + k')—^A^^-^ X log j . (52) 

We now turn to consider the contribution from taking in each wavefunction the terms that are linear 
in 5. These terms have been neglected in the literature. 

• Contribution proportional to 5\og{—Hr]) 

Again, as we are linearizing in 5, we take the (5-correction only from one wavefunction at the time. We 
start by considering the term in each wavefunction which contains the \og{—Hri): 

^''f' ^ i-, ^^TT7^4 (1 + e-'^" X 4 (-HV) ■ (53) 

Notice that this part of the Taylor expansion of (46), taken as isolated, does not converge for ry ^ 0, as 
the log(— i777) diverges. However we are now going to notice that the time integrals are dominated by 
rji ^ r)2 ^ Tjk where ry^ is the time of horizon crossing for the k mode {kr]}~ ~ —1), which implies that 
for 5 sufficiently small, wc can consider this isolated part of the expansion in (46). Once we perform 
the integral over time, the effect of this additional logarithm will be to multiply the expression in (48) 
by \og{—cH'qk), where c is some order one constant. This is so because the integrals in time in three 
dimensions read 

/O pr)2 
dm (ir/ir7?r/ie(*+*^i+*^2)^ie('=-*^i-^2)''% (54) 
-oo J — oo 

after performing the contour rotation, and it is straightforward to see that they are dominated by the 
times 771 ~ r/2 ~ 1/^ which is where the exponential suppression terminates. This is not changed when 
we add to the integral above a term of the form log(—i7 771^2)5 and therefore the result of the time 
integrations will be just to replace log(— if 771,2) with \og{—cHrjk), where c is some order one constant. At 
this point, the momentum integration is the same as before. In the limit (5 — > this term will contribute 
to the logarigthmic divergency by multipling the coefficient of order 1/5 of F{5) in eq. (49). This is the 
same term that in the former contribution was multiplying the factor log(A;/^). There are six factors 
of 5 \og{— Hr]kc) / 2 coming from the six wavefunctions integrated in the loops. Notice that there is also 
two analogous terms coming from Taylor expanding the factors of 0(771,2)'^ in the measure of integration 
around (5 = 0. In the (5 — limit, they also multiply the coefficient of order 1/5 oi F{5), with the only 
difference that this time they arc proportional to minus two factors of 5\og{—cHr]k)- Once we take the 
(5^0 limit, all of these eight terms give the same result as in (52) with the replacement: 

log(fe//x) ^ \og{-Hr^k) , (55) 

that is 

(C,Ki)C,^,(t))i-ioop,A,t^+oo,Mog = -(27r)¥^)(fc + fcO ^^^, ^^^3^ p X \og{-Hrik) • (56) 
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where we have neglected the divergent term and the order one constant c as it is irrelevant, and where 
the (5 log subscript represents the fact that this is the contribution from the terms linear in 6log{—Hrf). 

• Remaining (5— components 

We are now left to consider the contribution from taking in one of the wavefunctions the remaining 
term proportional to 5. Notice that this term is dimensionless and only a function of A;r/i,2 and ki^2'ni,2i 
depending on which wavefunction we are dealing with. So, if we were to perform the time integral, we 
would obtain an expression similar to (48). Then the remaining momentum integral, by dimensional 
analysis, would give a result proportional to 

Sk'^l^'^y Fi6), (57) 

where F(S) is dimensionless and divergent for S ^ 0. In order for a logarithmic running to come from 
this term, it would be necessary for F to have a double pole for 5 — 0. However, we saw in eq. (52) 
that the F{6) defined in eq. (49) had only a single pole, and it is quite straightforward to realize that the 
additional terms proportional to S that transform F into F do not increase the degree of divergency of F. 
In fact, it is quite tedious but straightforward to show that upon inclusion of the remaining (^-correction 
to the wavefunction in the time integrals, these are still dominated by 771 ~ ry2 ~ At this point, 

we look at the part of the (5-corrected wavefunction that we are considering, which now can be thought 
of as depending on the ratio ki_2/k (in the other case wc arc just left with a number of order one), and 
concentrate only on those terms that can induce a logarighmic running. These are only those ones that 
do not scale as a simple power law in ^1,2/^- The part we are considering of the (5-corrected wavefunction 
has a logarighmic dependence only in the IR limit k\^2/k — as \og{ki^2/k). However, in this limit the 
integral we started with in the limit of 5 ^ was power law convergent, and therefore the additional 
infrared log(fci^2/^) does not induce an additional logarithmic divergency. We conclude that F{5) has no 
double pole as (5 — >^ 0, and that therefore no additional logarithms appear when considering this part of 
the (5-correction to the wavefunction. 

Wc stress that these contributions proportional to 5 were forgotten in the literature [22, 21, 24] and 
led to parametrically incorrect results. 

Finally, we comment on the decision we took in eq. (44) to take the external wavefunctions directly 
in d = 3 dimensions. We can notice that if we had kept them in generic dimensions and we would 
have then expanded in 5 as we did for the other wavefunctions, then naively there would be terms 
that, multiplying Fq/5 from the loop integral would give rise to additional logarithm of the form \og{k) 
or log{H/fj,). However, if the divergencies are to be renormalized, than there is a counterterm that 
cancels the divergent term Fq/S of the overall integral. This counterterm multiplies the same external 
waveunctions as in eq. (44), canceling in this way all the additional logarithmic divergency that we would 
have had if we had kept the external wavefunctions in d dimensions ^. This justifies our rather intuitive 

^We will perform explicitly the renormalization in sec. 3.3. However, we can anticipate that schematically the 
counterterm will be of the form 

H~ ~ J (f+'x / a^+'{t) ^^^{d\f , (58) 
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approach of treating the wavefunctions that do not participate to the loop integral directly in three 
dimensions. 

Summarizing: in the calculation performed in dimensional rcgularization the effect of considering the 
terms proportional to 5 in the measure and in the wavefunction has achieved the result of transforming 
the logarithmic divergency log(fc/;u) found in the literature [22, 21, 24] into 

log{k/fi) + logi-Hrjk) - log{H/ii) , (60) 

where we have used that A; ~ — 1 . It is straightforward to see that the contributions we called B and 
C do not result in any logarithm and therefore we conclude that in dimensional regularization we also 
obtain: 

2 p2 tt8 i\j8 1 / H\ 

{cmkmi-ioo,,t-.^oo = -{2.f6('HUk')—^ . (61) 

3.2.1 \og{H/ij) without evaluating the loop integrals 

We would like here to show that, even without doing the actual calculation in dim. reg., it is possible to 
see that the simple fact that the integrals are dimensionally regularized forces the logarithmic running 
to be of the form \og{H/ In order to this, it is however necessary to assume that the ^ two-point 
function does not depend on time. It is possible to show this without actually having to do any loop 
integral, as we will show in sec. 5, or, by actually doing the calculation in dim. reg. at finite time as we 
will do in sec. 3.2.2. For the time being, we will assume this to be the case and we will therefore take the 
extremum of integration of the time integrals in the loops to plus infinity. 

At this point, by proper redefinition of the variables of integration, it is possible to make the loop 
integrals dimensionless functions of 5 only, and from this extract the form of the logarithmic running. In 
fact we can take for example the loop term in (42) that we reproduce here for convenience: 

(C,Kt)C,^,(i))i-ioop,A = -(27r)3<5(3)(fc + P) I28ci^^i^p (62) 

X T^kim) T^k'im) <(m)<*(^2) <(m)<*(^2) , 

where the number of derivatives comes from imposing that the counterterm has the same powers of the cutoff Ay 

as the terms it renormalizes, as it is standard in non-renormalizable effective field theories, and where we have 
not been careful in distinguishing spatial or time derivatives, as in general all will appear in the counterterms. By 
treating (58) perturbatively in (14), we will obtain terms schematically of the form 

The first two wavefunctions on the left are the same external wavefunctions we get in the loop. Taylor expansion 
of the integrand in 5 ^ (and using eq. (46) since the integral is dominated by krj' ^ —1), shows that no addi- 
tional logarithmically divergent terms linear in 6 are generated by the integrand. This implies that the additional 
logarithmically divergent pieces that comes from considering the external wavefunctions and the i5-function of the 
external momenta in 3 -|- 5 dimensions get canceled by the counterterm. 
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where we have passed to conformal time, taken the rj ^ hmit and we taken the external wavefunctions 
in three dimensions. We can make all the quantities within the integral dimensionless by multiplying or 
dividing by the external momentum k. This is possible because the scale factor is unchanged in dim. reg: 

and because the wavefunctions in dim. reg. are just functions of kr] apart for multiplicative factors: 

^fei (^l) = ^ , (3+S)/2 3+5)/2(-^l^l) (64) 

2 (-HM^y^ ' k\ 

( (h\ \(3+5)/2 
2 fc(3+5)/2(^^/^)(3+5)/2 ^(3+5)/2 1^ ) ^ 

~ / . \l/2 Z,(3+5)/2"-^ I U ''^''1 i ' 

2 [-HMl.ix^) / >f \k J 

where G{k\/k,kir]) is a dimensionless function. By extracting the relevant powers of k and H, the 
remaining integral is just a dimensionless function of 5. We obtain: 

(4W4-,W>:-,.P,A = -P-)'*<^'(* + *')128<|^-5^pX_^/W (65) 

5 



where I{5) is a dimensionless function of 6. In the limit (5^0, the factor of {H/jjiY can be Taylor 
expanded to give 



This quantity will multiply any pole coming from the integral I{8) as (5 — 0, giving rise to a logarithm 
of the form 

This shows that i/ there is a logarithm as a result of this contribution then it is of the form log(iJ//i). It 
is straightforward to sec that the same result holds also for the other diagrams. This is an alternative 
proof of the fact that the running has the form \og{H/ jjb) that does not rely on us being able to actually 
perform the loop integrals. 

3.2.2 Calculation in dimensional reguleirization at finite time 

It is not difficult to generalize the calculation we did in sec. 3.2 to finite external time. Apart from the 
fact that the integrals becomes slightly more complicated the only subtlety is in how to deal with the 
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corrections from the measure and the wavefunction that are proportional to 5 (see eq. (46)). In sec. 3.2 
we saw that, when the external time is taken to be plus infinity, the time integrals were dominated by 
tji ^ r]2 ^ rjk, where r/^ is the horizon crossing time of the mode k. This allowed us to replace the term 
5log{—Hr]i^2) in the wavefuntions and in the measure with log( — Hr]k) and argue that the remaining 
(^-corrections to the wavefunctions did not induce addional logarithmic divergencies. It is clear that if 
the external time r) is late enough so that the external mode k is outside of the horizon (i.e. —krj <^ 1), 
this still applies. In this regime the (5-corrections can be treated exactly as in the case of infinite external 
time. 

A straightforward calculation gives: 

(C,-(i)C,-,(0)i-ioop, = -(2.)3^(3)(fc + fc-0^^i^pxlog(^-j (68) 



X Re 



-|(z + krjf (4 + kr] {2i + kri) (4 + kri{6 + ikr]))) 



As we will discuss in detail later in sec. 5 the loop integrals performed at finite external time are, at 
least in principle, more regular than the ones done at infinite external time. This is so because if we keep 
the external time finite the most ultraviolet k modes are inside the horizon and therefore the integrals 
receive additional oscillations in the UV that improve their convergence. It is therefore possible that the 
logarithmic running that we found in the infinite external time limit disappears in the calculation done 
at finite time. It would be replaced by a time-dependent term that diverges as log(?7) as we send the 
external time r] to zero. We find that this is not the case. We will further explicitly verify this in sec. 3.3 
where we will find the counterterms that allow us to reabsorb the divergencies. This proves that the ( 
correlation function is time-independent at one-loop, and it also shows that the logarithmic running we 
computed in (68) can indeed be computed without need of renormalizing the theory. 



3.3 Renormalization 

We now proceed to the actual renormalization of the correlation function of C(*)^ for the case of the 
large tt^ self-interactions. This theory is so simple, that this process is not complicated. Notice that this 
interaction is irrelevant, or equivalcntly non-renormalizable. 

In non-renormalizable effective field theories, the divergencies are meant to be absorbed order by order 
in an expansion of energy over cutoff, or equivalently derivatives over cutoff. The effect we computed in 
(73) is from a dimension six operator (in canonical normalization), and therefore it is suppressed with 
respect to the tree-level result by (H/Au)'^, where Af/ is the unitarity bound defined in eq. (40). If we 
work in dim. reg. only the logarithmic divergencies appear and therefore we have to look for dimension 
eight quadratic operators. There are three dimension eight quadratic operators compatible with the tt 
shift symmetry [3]. The counterterm Lagrangian is therefore given by: 



Sc = J d'^x a(tf+^HMl^ 



M8 



(69) 



where the Ci, C2, C3 are dimensionless numbers, M^/{HMp^)^ represents the cutoff suppression, and 
we have extracted a factor of HMp^ that represents the normalization term of the two derivative kinetic 
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term ^. 

Notice that the procedure would be different if we instead decided to perform the renormahzation 
using the cutoff regularization. In this case, all the lower dimension operators would be generated with 
some power law divergency. These include the lowest dimensional one: the wavefunction renormahzation. 
Notice however that the coefficient of the term (diir)'^ is fixed by symmetries to be HMp^, due to the 
fact that TT is the Goldstone boson of time translations [3]. It cannot be renormalized. This is quite an 
unusual non-renormalization theorem that affects this term. This means that the power law divergent 
terms will only renormalizc the tt^ term out of the two lowest dimension one. This effectively induces 
a speed of sound of the fiuctuations Cg that is different from one, unless it is tuned to be equal to one 
(and this is technically possible since we are dealing with a power law divergency) [3, 14]. This was 
expected, as typical in effective field theories all allowed operators tend to be generated suppressed by 
the same scale unless we perform some tuning or unless there is a symmetry protecting it. In this 
case, there appears to be no symmetry protecting Cg = 1 and we decide to tune it. Notice that in 
this paper we are concentrating on the conceptual issues arising from the study of loop corrections to 
inflationary observable, and we consciously decide to study a consistent, though tuned, theory to simplify 
the calculations. 

Connected to this, we notice that there are higher derivative quadratic operators of the form {d^7r)^/Afj 
which would naturally give a correction to the ^ two-point function larger then the one in (73) by a factor 
of the order (Au/H)^ » 1. Again, we are assuming that these terms are tuned away. 

Coming back to reabsorbing the divergency of the dim. reg. loop result of (73), the counterterms 
give a correction of the form 



(C^(*)) counter term — —2H^lm 



J—OO- 



(70) 



(2) (2) (2) 

where H-^j. = —Lc , where Lc is the counter term Lagrangian. This gives: 



1 1 

(a(t)a'(i))cour.ter term = (2^')'^^'^ + ^ ^ (71) 

X [- (Ci + 3C2 + C3) (1 + (krif) - 4 (3Ci + C3) (kr,)^ + 2{C\ + C2 + C^){kiif] . 

These terms have to cancel the pole in \/8 coming in the loop integral. At finite time, this is given 
by expression (68) with log(if/ ji) taken to be equal to one. Notice that the finite time results depends 
on (krf)'^, {krf)^, {kr]y,{kri)^: imposing the cancellation of the poles leads to four equation in three 
unknowns: Ci,2,3- The solution still exists: 

Ci = 47rci^, C2=^^ncf-, C3 = -yTrc^l , (72) 

providing another non trivial check of our results. The pole 1/6 represents the fact that these counterterm 
are actually divergent in three dimensions. To each one of the coefficients we can add a finite term that 



''In the unitary gauge where the tt Lagrangian is usually constructed [3], these term arise for example by 
operators such as {d^Sg^'^)^, {dtdiSg^'^)^ , (dtSKl)^, where Kl is related to the trace of the extrinsic curvature of 
uniform inflaton surfaces (see [3] for details.) 
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can be determined by imposing the two point function at some renormalization scale to be equal to some 
quantity. We do not perform this last step here, as it is not the main point of our work. 

In fact, our result shows that the divergency we found in sec. 3.2 are reabsorbable by local couterterms. 
This is yet a third independent proof of this (the first one was obtained by performing the momentum 
integration first in the loop integral regularized with a cutoff; the second one was by performing the dim. 
reg. loop integrals at finite external time in sec. 3.2.2; we will give a fourth proof of this in sec. 5). 

Though the finite terms in (72) could be determined only by performing the loop integrals in sec. 3 
paying attention to the finite terms it is worth mentioning that as usual the effect we have computed can 
be much larger than the finite effects if we decided to renormalize the theory at an energy scale fx very 
far from H. This is actually what makes sensible to compute the running and not the finite terms. 



3.4 Summary of calculation with large tt^ self-interaction 

The case where the inflaton has a large self-interaction of the form tt^ has allowed us to calculate one- 
loop corrections to the inflaton two-point function. We have performed the calculation in two different 
kind of regularizations: one with a sharp cutoff in momentum and frequency space, and the other in 
dimensional regularization. We have also performed the dim. reg. calculation at finite time, and we have 
also shown the form of the logarithmic running by simply making dimensionless the loop integrals. We 
have then performed the renormalization of the correlation function, by explicitly finding and solving for 
the relevant counterterms. 

We have found that all the calculations agree. The result is: 

(C,-(OC,-,(t))i-ioop,t^+oo = -{2.f6('Hk + k')^,j^y, X log ^- j . (73) 

As we had anticipated the logarithm is \og{H/ We find this to be a very sensible result. First of all, 
contrary to the log(fc//x) found in the former literature, the real space version of (73) is symmetric under 
the rescaling 

a— ^Aa, x^x/\, k^Xk, (74) 

a symmetry of the problem. Second, it makes sense from a physical point of view: the Hubble scale 
is cutting off the infrared behavior that would otherwise be there in Minkowky space. The resulting 
logarithm is of the form log{H/fi), which is similar to the form found in scattering amplitudes: log(£'/;u) 
where E is the center of mass energy and /j, is the renormalization scale. The energy probed by the 
interactions during inflation is of order H. 



4 Gravitational interactions with massless scalar fields 

We now turn to the inflationary theory that Weinberg originally studied in [2]. This is a theory where 
an inflaton with a standard kinetic term is rolling down a flat potential and is interacting gravitationally 
with TV" massless scalar fields. The Lagrangian is of the form: 
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n=l 



(75) 
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Though we are considering purely gravitational interactions, the contribution from the N massless scalar 
fields running in the loops will be enhanced by a factor of N with respect to the analogous interac- 
tions coming from the graviton and the inflaton running in the loop. For this reason we can avoid the 
complication of letting those run in the loops and concentrate on the a scalar fields ^. 

We do not need to redo all the calculation. It is straightforward to include in the results of [2, 21] 
the correction coming from the terms proportional to S. In fact the only relevant terms proportional to S 
were those of the form log{—Hri) coming cither from the wavefunction or from the measure of integration. 
After performing the time integrals, which converge and are dominated by r]i^2 ~ Vki those terms become 
of the form log(— -ffr/fc). The contributions with two insertions of the cubing interaction Hamiltonian are 
of the form: ^ 

k"" (^^) F{S) X (1 + Slog{-H7^,)) , (76) 

where D is the correct number of dimension for the integral in three dimensions, while the overall factor 
of comes from the two measures of integration times a (3 + 5) dimensional 5-function. As before, F{6) 
will be a dimensionless quantity of the form 

Fid) = ^ + F,, (77) 

for S ^ 0. The coefficient of the term Slog{—Hr]k) is equal to one because this contributions involves six 
wavefunctions and two measures of integration. Eq. (76) can be expanded for small S to give: 

k' (l + 6log +dlog{-Hr,kij + Fi) ^ k' (l + dlog ) (^^ + Fi) , (78) 

which means that the coefficient of the logarithmic divergency is of the form 

log(^) . (79) 

This analysis can be easily extended to the diagrams involving an insertion of the quartic interaction 
Hamiltonian. In this case the integral in time converges and is dominated by the time of the horizon 
crossing for the k-mode. The result after the momentum integration will be of the form: 

A:^ (^)'^('^) X (i + Slogi-Hrik)) , (80) 

which is of the same form as (76). The factor in front of S\og( — Hf]k) is equal to one and this time it 
comes from the fact that there are four wavefunctions contributing as 26log{—Hrik) and one measure of 
integration contributing as —Slog{—Hr]k). [21] shows that in the case of N spectator scalar fields, F{S) 
does not have a pole at (5 — ^ 0, and therefore no logarithm can come from this term. 

^The calculation done originally by Weinberg has a rather minor numerical mistake due to an inconsistent 

implementation of the ie prescription. This was noticed and fixed in [21]. Both [2, 21] performed their calculations 
in dimensional regularization and did not include the contributions proportional to S from the scale factors in the 
measure of integration and from the rf-dimensional wavefunctions and therefore obtained an incorrect logarithm of 
the form log(fc//x). 
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Summarizing, in the case of A^" spectator scalar fields the form of the logarithm is \og{H/jjL). By 
applying our correction to the results of [2, 21] we find: 



(81) 



4.0.1 Tadpole diagrams and non IPI diagrams 



In computing one loop corrections to correlation functions of C, there are important tadpole diagrams that 
have been so far neglected. These tadpole diagrams arc important not only for the one-point function, but 
also, by attaching them to propagators, for all correlation functions. Because of translation invariance, 
the external line attached to a tadpole diagram has to have zero wavenumber. In the former example 
(tt^ interaction) these diagrams were therefore zero because the field tt had an exact shift symmetry so 
that every tt had a derivative acting on it. For the spectator scalar fields this is not case. The tadpole 
diagrams have to be included. The expectation value of Q in this theory is given by: 



(Cfe(i))i-ioop = -2Im 



(82) 



= -(27r)=^(5(3)(jt)Im 
where we defined 



£ dt,a{hnt{t) Ut{ti)\ (2(p.(ti))o - 3(p.(ii))o) + ^(p.(ti))c 



{Pa{t))o = { 1 , 



(83) 
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and where the expectation value cannot depend on x because of translation invariance. At this order, 
Hint = -4nt and is given by: 
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We have assumed that the expressions for (/Oo-(t))o and {pa{t))o, which are UV divergent, have been 
renormalized and have been made finite (at most dependent on the renormalization scale). Their actual 
value is not important. The wavefunctions of C and a in three dimensions are: 
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1/2 ^,3/2 
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(85) 



A diagrammatic representation of this calculation is given in Fig. 1, where the distinction between 
dashed and continuous lines will be explained in sec. 5.1. The terms inside the Imaginary part in (82) 
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da 

X ;'' X 

da 

Figure 1: Tadpole diagram. Dashed lines represent correlation functions, continuos lines represents 
Green's functions. We will explain the origin of this notation later in sec. 5.1. 



are not real and therefore the diagram is not zero. However, the wavefunction in (85) diverges for k = 0, 
which means that the tadpole diagram is actually infinite even after having renormalized (pcr(^))o and 
{Pait))o. This is a problem that arises because eq. (13) does not fix the gauge at A; = 0. For tree-level 
calculations this is not important, as the zero mode of <^ is not observable. However, the tadpole diagram 
can be attached to a (^-propagator and affect at one-loop also correlation functions with k ^ 0. Notice 
that the trilinear Lagrangian in (, Ci^s that can be found in [1] does contain vertexes of the form C(^C)^ 
and C^C which allow for a A; = mode to be attached to a propagator and give a non-zero (and actually 
infinite) result. 

However this problem can be dealt with in a rather straightforward way by ensuring that the zero 
mode of C, is zero. This can be done by realizing that the tadpole is proportional to the expectation 
value of the stress-tensor of the a fields in the unperturbed metric. Since translation invariance forces 
the expectation value to depend only on t, this expectation value does nothing other than changing the 
background solution. The correct way to deal with the C, zero-mode is therefore to define as in equation 
(13), but where in this case the background quantities (and in particular the scale factor) satisfy the 
following equations: 

^M^,H^ = ^^{tf + V (<^(0) + {Pait))o , (86) 
M|i {3H^ + 2H) = -^^{tf + V {4>{t)) - {p,{t))o . 
Formally, this corresponds to manipulating the Lagrangian of (75) in the following way: 




(Nf Nf 
1=1 i=l 

The background solution comes from imposing that the first line of (87) starts quadratically in C (which 
just means solving the unperturbed equation of motions coming from the first line). Instead, the second 
term in the second line ensures that the contribution coming from contracting the two cr's in the first term 
in the second line (i.e. every tadpole diagram or every subdiagram containing a tadpole subdiagram) is 
cancelled. This is represented pictorially in Fig. 2 for the one-point and the two-point function. 
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da 




=0 



-XX- 



-XX- 



Figure 2: Upper: Cancellation of the tadpole diagram. Lower: Cancellation of the tadpole subdiagram 
in a two-point function diagram. Dashed lines represent correlation functions, continuos lines represents 
Green's functions. We will explain the origin of this notation later in sec. 5.1 



5 Time independence of ( out of the horizon 

In the former sections we showed that the one-loop corrections to the two point function of give rise 
only to logarithmic running of the form \og{H/ ji), where (i is the renormalization scale. As part of the 
calculation we showed that loop corrections do not give rise to any additional time dependence of the 
two-point function of C,^, which therefore stays constant once the mode k is outside of the horizon. There 
is an important subtlety that was stressed by Weinberg in his original work [2, 22]. It is associated with 
the fact that in order to compute the logarithmic running we did not need to renormalize explicitly the 
theory. In other words, we had simply to assume that all the divergences could be reabsorbed by diff. 
invariant local counter terms in the Lagrangian. It is possible that this might not be the case for the 
following reason. If we take the calculation we did in dimensional rcgularization in sec. 3.2 we notice that 
first we performed the time-integrals in the loops by sending the external time to plus infinity and then 
second we performed the momentum integral finding a divergency that led to a logarithmic running. It is 
possible that this divergency might not be absorbable by a diff. invariant counter term in the Lagrangian. 
This could happen if in doing the time integrals instead of sending the external time to infinity we were 
to keep it finite and the divergency only appeared as we sent this time to infinity. We are unable to say 
anything about this unless we perform explicitly the renormalization. 

We have not ignored this issue completely in the former sections. For example when we did the 
rcgularization with a cut-off we first performed the momentum integrals at finite external time and 
second we performed the time integrals and in doing this we were able to take the external time to 
infinity only at the very end. Since we found the same result as in dim. reg. this proved that the C, 
correlation function could not depend on time and that the divergency we found should be re-absorbable 
by a counter term. 

The only issue with this second calculation is that the cutoff in frequency and momentum space is 
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not diff. invariant, which might lead us to suspect either that the logarithmic divergency that we found 
is not correct, or that there are additional time-dependent factors that were lost because of this. We 
believe that this is not the case, as the logarithmic running is an infrared quantity that should not be 
sensitive to the details of the UV regularization, particularly in our case. 

In order to further address this point, in sec. 3.2.2 we performed the calculation in the dim. reg. for 
finite external time. We found that the divergency was unaltered. This showed that the divergency was 
real and had to be rc-absorbed by a counter-term. 

The calculation at finite external time we just mentioned was done only for the vr^ theory. It is 
pretty clear that it can be extended to the case of the spectator a fields, but we have not done this 
explicitly. Here we will offer a proof of the constancy of C, out of the horizon directly for the theory 
with N spectator a fields. We will then generalize it to the theory with the large vr^ self-interaction. 
We find that this additional proof improves our intuition on loop corrections to infiationary correlation 
and it will further guide us in speculating on how our results are generalized to higher loops and to 
different interactions, which we will do in sec. 6. This way of doing the calculation can be thought of 
as doing the loop calculation by first doing the momentum integrals and then doing the time integrals. 
Since if one does the calculation in this order there are no subtleties with the possible appearance of fake 
UV divergencies we can assume that the quantities we will deal with after the momentum integral are 
directly the renormalized, physical ones. This is why this method of solving the loop will be so intuitive: 
it will reduce to solving the linear equations of C, in Fourier space working directly with the renormalized 
stress tensor. Though very intuitive this method will be not be powerful enough to show that there is 
logarithmic running in the two point function and calculate its coefficient. But it will make very clear 
why the C, correlation function is constant in time and scale invariant once the mode is outside of the 
horizon. 

As wc did in the introduction, we would like to stress that it would be a tragedy for infiation if the 
correlation function of or of some other related operator did not become constant in time after horizon 
crossing because of some loop interactions involving modes of Hubble size. If this were to be the case, 
than we could expect the same to happen for interactions during the unknown part of the history of the 
Universe from the end of infiation to nucleosynthesis during which the modes relevant for observation 
are well outside of the horizon. In fact, as it will become clear later in this section, the quantum nature 
of the loops is not really relevant here: any type of fluctuations can lead to the same effect, and it is 
therefore extremely important to show that fluctuations on scales smaller or comparable to the horizon 
cannot influence modes outside the horizon. It is the constancy of C during the unknown parts of the 
history of the Universe that makes inflation, or any other theory of the early Universe, predictive. 

Summary of the section. In order to prove the constancy of C, we will introduce a alternative way 
of looking at loop corrections to correlation functions. This will take some time and therefore we would 
like to anticipate here the main logic and results. Concentrating on the theory with N massless spectator 
scalar flelds a, let us start by noticing that the operator at time t is given by an expression of the 
following form: 

Ut) = [ dt' Gf^\k, t, t')T,,<:;^,{k, t') . (88) 

J —oo 
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Here T(^a)iJ.u is the stress tensor of the a fields: 



_ (89) 

and the subscript represent the fact that in some cases the operators a are to be computed in the 
presence of a fluctuations. G^'^^ {k, t, t') is the retarded Green's function for C, associated to the T(^a)fiu- 
Since the a's interact with the inflaton only through gravitational interactions, it is Ti^a)nu that sources (. 
The above expression just comes form inverting the linearized Einstein Tensor (which is invertible once 
wc fix the gauge), and in general it will be a complicated, apparently non-local expression. We will give 
the explicit form for later in this section. Since we are interest in proving that correlation functions 
of (^k become constant when the mode k is outside of the horizon, we will be interested in considering t 
as a late time after horizon crossing as possible, for example the end of inflation. If t* is some time after 
which the mode k is outside of the horizon such that: 

^ = eout < 1 , (90) 



a{t*)H 
we can rewrite (88) as: 

Ck{t) = f dt' Of'^'ik, t, t')T,,cMk, t') + f dt' Gf-'^^ik, t, t')T,,c-Mk, t') ; (91) 

J-oo Jt* 

The first term on the right-hand side picks up the contribution to Cfc(^) generated by the sources that 
acted up to the time t*, i.e. while the k mode is inside the horizon and as it goes out of the horizon, 
while the second term represents the contribution to Cfc when k is already outside of the horizon. We 
will analyze each of the two terms separately and we will show that they contribute to the C two point 
function in a way that is scale invariant and time independent. Let us anticipate the main result here. 

The contribution from the first term on the right-hand side can be thougth of as the free evolution 
of C as determined by some initial condition at time t* . Because the free evolution of the Cit makes it 
freeze after horizon crossing this means that this contribution will become time-independent after horizon 
crossing. We will argue that this implies that the contribution is also scale invariant ^. 

The contribution from the second term on the right-hand side of (91) is instead a bit more complicated 
to deal with. It represents how the fluctuations in the T(^a)fiu source the Cfc mode when this is well outside 
of the horizon. Since the u's have no potential term, their fluctuations are relevant only for the a modes 
smaller then the horizon. Still, the superposition of two short scale modes can lead to a long scale mode, 
and therefore affect Cfe- 

Concentrating on the two point function, there are two kind of contributions from these terms. 
The flrst, that we will call cut-in-the-middle diagrams, represents how the free two-point function T(^a)tiu 
affects the (k two-point function. We will see that once the mode k is well outside of the horizon, the 
various Hubble patches spanned by one wavelength of Cfc become uncorrelated. This means that the T(^a)nu 
fluctuations quickly average out and are not able to provide the coherent effect that would be necessary 
to source on large scales. Cfe becomes therefore time-independent once outside of the horizon, and 



^We will work initially in the approximation of de-Sitter space, and we will then generalize to deviations from 
exact de-Sitter. 
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again the symmetry under rescaling of the scale factor a will force the contribution from these terms to 
be scale invariant. 

The second kind of contribution from the right-hand side of (91) comes from computing the expec- 
tation value of T^a)^iu the presence of a background (k- The background (^k can be correlated with a 
freely evolving (^k to obtain a two-point function. We will refer to this contribution as the cut-in-the-side 
diagrams. The reason why this contribution shuts down as the k modes go outside of the horizon is that 
the perturbation to the T(^„^^y due to the presence of a background C, goes to zero. This happens because 
the background Q at this order in perturbation theory evolves freely. When a free C, mode is outside of the 
horizon it becomes constant and a simple redefinition of the scale factor a which is locally unobservable. 
This means that the C, mode cannot induce any physical perturbation to T(^a)iJ.u- This implies that the 
effect of this term is again time independent, and due to the usual symmetry under rescaling of a it is 
scale invariant. 

This will conclude the altenative proof that ^ is time-independent and scale invariant out of the 
horizon for the case of the N spectators a fields. Finally, we will be able to comment on slow roll 
corrections to our calculation and extend this particular proof to the case of the large tt^ self-interaction. 

5.1 An alternative diagrammatic expansion for loop corrections 

We start the second proof of the constancy of when it is outside of the horizon by presenting an 
alternative way of organizing the calculation. This approach was originally developed in [29] for a 
restricted set of theories, and it was noted in [21] that the derivation was not consistent with the ie 
prescription for choosing the interacting vacuum in the past. Here we will generalize the approach of [29] 
to more generic theories and we will show how the correct i e prescription can be implemented. 

For concreteness let us specialize to the ( two-point function. We start by taking expression (14) and 
inserting the unit operator 

Uintit, -f)Url (t, -f) , Uint{t, -f) = Te-' ^-t ^^"*(*') , (92) 

between the two C's, to obtain 

(C'W) = {{Ur^tit^ -oo-)at)Uint{t, -f)) {Urlit, -f)C,{t)Uint{t, -oo+))) , (93) 

Here T is some arbitrary early time in the past, and we stress that no rotation of the contour of integration 
is performed for this operator. Ignoring for a moment the issue of the i e prescription, we have written 
the expectation of the operator C(i)^ as the product of the two freely evolved C(0's each evolved with the 
interaction picture time evolution operator Uint- In other words, the C(i)^ correlation function is simply 
given by the correlation function of the evolved C(^)'s- A closer look at (93) might let us think that this 
picture does not quite work because of the difference in the extremes and the contour of integration of 
the time integrals in the various evolutors. In order to make it clear that these differences do not play 
any significant role, we can deform the contour of integration for the two external Uint to obtain: 

{C\t)) = {U^nli-T, -oo_) {Url{t, -f)at)Uint{t, -f)) {Url{t, -f)C{t)UUt, -f)) UU-f, -oo+)) . 

(94) 

The two most external Uint have the function of projecting the free vacuum into the interacting vacuum, 
while the rest of the i/'s evolve each from the time —T to the time t. Notice that, in this approach. 
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we cannot take the time —T to — oo from the beginning of the calculation because we would lose the 
capability of projecting on to the true vacuum. We will comment shortly of the physical meaning of this 
point. 

Going back to eq. (93), we can Taylor expand in Hint to obtain 




dtN- 



t2 



dU 



/tjV rt2 
dtN-1 ■ ■ ■ dti 



Hintih), Hint{t2), ■ ■ ■ [Hint{tN),C,{t)\ 
Hintit'i), Hint{t2), ■ ■ ■ [Hintit'^), C(0] - 



-OO- 

-f 



-OO- 

■T 



-OO- 



-T 



J -T 



-T 



where ^ ^ ^ 

J dt' [Hint{t'),C,{t)\ lf~ = / dt' Hint{t')m - t _dt' at)Hint{t') 

Expanding (95) up to second order in Hint-, we obtain 

{eit)) = {e{t))cis+{e{t))ciM , 

where we have defined 

{C\t))cis = -2 Re 



(96) 
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(98) 



The subscript ciS denotes what we call cut-in-thc-sidc diagram, while ciM denotes cut-in-the-middle 
diagram. If we remind ourselves that the ( retarded Green's function is given by 



Gfix,x') = i0it-t')[ax),ax')] , 



(99) 



where the commutator on the right is taken on the free fields, we can identify the structure we anticipated 
in eq. (88), with the additional subtlety that each term of the commutator is taken on a different time 
path (see eq. (98)). Alternatively, as in eq. (93), one can deform the contour of integration for the integral 
that goes up to — oo, in such a way that it goes to — T on the real axis, and then goes from there to — oo 
on the imaginary axis. In this way, we can treat the integration from — oo to —T on the imaginary axis as 
the projection of the state at time —T onto the interacting vacuum, and we can interpret the diagrams 
in (98) as representing the evolution of the operators from time —T to time t starting in the interacting 
vacuum. This subtlety will be largely irrelevant in proving the theorem. 

By specializing to the case of the N spectator a fields, we can see that there are two ways in which the 
equation solution for ( can be perturbed by interactions of the form (a^. The first, which corresponds 
to the CIM diagrams, is by considering vacuum fiuctuations of two a fields, that combined in a T(^(j),iJ.v 
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Figure 3: Cut-in-the-middle diagrams. Continuos lines represents Green's functions, dashed lines repre- 
sent free fields, and crosses represent correlations of free fields. Two crosses have to be contracted together 
in order for the diagram not to be zero. 



generate a C mode which then propagates freely. The correlation of two of these terms is the CIM 
diagram represented in Fig. 3. The notation is such that a dashed line corresponds to a free field, while 
a continuous line represents a retarded Green's function. The crosses represent correlation of free fields 
(two crosses have to be contracted together in order for a diagram not to be zero). 

The CIS diagram represents instead the correction to the C, correlation function due to the pertur- 
bation to the evolution of a primordial C, fluctuation caused by the interaction with the a fields. It is 
represented in Fig. 4. We see that an original C, vacuum fluctuation interacts via the ^cr^ interaction with 
a cr vacuum fluctuation and modifies that fluctuation. This fluctuation evolves for some time and then 
interacts again with another a vacuum fluctuation, generating a Q fluctuation that then propagates freely 
up to time t. The correlation among the original a vacuum fluctations and of the original C, fluctuation 
with another C, fluctuation that instead propagates undisturbed up to present gives rise to the diagram. 
Notice that the subdiagram on the left of the last C, Green's function is nothing but the expectation value 
of the T(^a)fii/ time ti in the presence of a background (" mode. 

The names cut-in-the-side diagram and cut-in-the-middle diagram are meant to point out the fact 
that in the CIM diagrams both of the ^ operators undergo a perturbed evolution, while in the CIS 
diagram it is only one of the two C's that is perturbed. Though here we have concentrated only on 
the IPX diagrams with insertion of the cubic interaction Hamiltonian, it is straightforward to include 
also the quartic and the non IPI diagrams. Clearly, the CIM and the CIS diagrams are physically 
well distinct: for example, the CIS diagram is proportional to the primordial ( vacuum fluctuations, 
and goes to zero in the limit that there are no C primordial fluctuations. On the other hand, the CIM 
diagrams are independent of the primordial ^ fluctuation. In generic situations, these two diagrams can 
give parametrically distinct results. However in the case we are considering here where the a and the 
( flelds are massless and undergo quantum vacuum fluctuations the two terms are comparable. In fact 
in the interacting theory the a and ( fields are mixed and both diagrams give equal size contributions. 
Equivalently although the overall result is independent of the time T that dependence only cancels in 
the sum of the CIS and the CIM diagrams. 

We will find the organization scheme and the physical intuition provided by these diagrams extremely 
useful for proving the time-independence of the Cfc correlation function once the mode k goes very well 
outside of the horizon. Furthermore in sec. 6 the intuition we will develop will allow us to generalize our 
results in various ways. 
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Figure 4: Cut-in-the-side diagrams. Continuos lines represents Green's functions, dashed lines represent 
free fields, and crosses represent correlations of free fields. Two crosses have to be contracted together in 
order for the diagram not to be zero. 



5.2 Early time contribution 

Let us now begin to analyze the two contributions from (91), and let us start from the first term on the 
right-hand side, that we can rewrite as 

CkAt) = f dt' Gf''"'{k,t,t')T^,)^,^,{k,t') , (100) 

where the subscript i has been introduced to represent that we are talking of the contribution to ( from 
the first term on the right-hand side of eq. (91). Notice that since in this subsection we will automatically 
consider the sum of the CIM and CIS diagram, the subtlety about the two different time-integration 
paths in the Green's function in (98) is irrelevant. In (100), the source is allowed to act only for the 
time up to t* when the mode k begins to be well outside of the horizon, and therefore represents how 
is affected by early interactions. But exactly because of this it is easy to realize that as far as this 
contribution is concerned the source will just create a certain Ck{t*)j Ck{t*) at time t* which will serve as 
initial conditions for the subsequent free evolution. Therefore, this term is equivalent to study the time 
dependence of a freely evolved ({t) given some initial condition Cfe(**)) Cfc(**)- The solutions in exact 
de-Sitter space are 

a(0 = CM(0«it + c£*(*)4' (101) 

^°We will generalize later our results to deviations from de-Sitter. 
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where the wavefunction (^f is given by: 

a[,{t) = A,{C\C,k)^{kvcos{kr,)-sm{kri)) + A2{C\C\k)^{^^^^ (102) 

where ^i(C*5 C*; k) , and A2{C*, C*) k) are the two integration constants determined by and at time 
t*, and 77 is conformal time. Independently of the value of ^1,2, the solution for (^^ goes to a constant 
value exponentially fast in cosmic time t. 

The symmetries of de-Sitter space suggest that we demand that the contribution to the variance of 
^ in physical real space that is given by modes between two fixed physical scales (both smaller that 
H~^) should be independent of time. Satisfying this constraint implies that the coefficients Ai^2 should 
be independent of k which then results in a scale invariant spectrum for (. We thus expect that the 
k dependence of Ai^2 is slow-roll suppressed and that the power spectrum of ( is very close to scale 
invariant. 



5.3 Late time contribution: cut-in- the-middle diagrams 

Let us now analyze the contribution from the second part of eq. (91), that we denote with the subscript 2' 

Ck,2{t) = f dt' Gf^\K t, t')T^,)c,,,u{k, t') . (103) 
J t* 

The contribution from these terms represents how the small scale tr's affects (" even when the mode k is 
well outside of the horizon. It is worth stressing that T((^)^(/c) is not zero in this regime. Two small scale 
a modes can combine to create a low-fe T((^)^(fe). 

Let us start with the CIM diagrams, where T(^^^ is evaluated in the unperturbed metric. In this 
case it is straightforward to find the equation of motion relating the a's to ( because the cr's live in 
an unperturbed metric. Further, since, as it will become evident soon, the CIM diagrams involve the 
self-correlation of four a's, we can neglect all the terms of order O (C^"^) from the equations of motion. 
The simplest way to find this equation is to use the action we wrote earlier in (84) at order O (C^^^) and 
to derive the equation of motion for (. Alternatively, one can solve the (0, 0) and (0, i) Einstein equations 
to find N and A'"' at order ^ and (but not at order (^a'^), and plug these back into the continuity 
equation 

V^T'*'° = . (104) 

One finds: 

and the resulting equation for is 

h\ ■ 92 1 / ^ {dia 



(+(^3/r-2- + jJC-^C-j^(^.*^ + ^l (106) 
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Notice that in all of the above expressions (105) and (106) the terms of 0{a'^) should be understood as 
having the zero mode subtracted, a result of the tadpole subtraction procedure. 
The above equation can also be written in terms of T(^c)iJi,v'- 

H h\ , d'^ ^ 1„ lfdt3^^HH\fl^^ , ^ 



(107) 



We are interested in solving this equation when the mode A; of ^ is well outside the horizon. We can 
therefore expand in powers of the external momentum (keeping track of the difference between di{ad'^a) 
and di&d^a), and keep only the leading terms. We obtain: 

Notice that at late time there are two sources for T^^, which represents the energy density of the 
(j's, and ^diT^\ which represents a non-local term proportional to the divergency of the velocity of the 
(j's. 

Eq. (108) can be integrated to give: 

where we have used the fact that we are interested in the contribution of the source from onwards. 
Notice how this expression is indeed of the form of eq. (103). 

We are now ready to begin the computation of the CIM diagram. We have: 

/ ^/'+„^3- 



{c.,mMt))c,. = I -ih I (^) (i - ^Ij^) (i - X (110) 



'PI 

There are six correlation functions to analyze. Let us start with the one involving 

{{a\m\it')) , (111) 

We notice that this operator is divergent, and here and in the rest of the section, the operators appearing 
in the correlation functions should be understood to be regularized. Notice that here we are not interested 
on the numerical value of these operators, but just on how they affect the late time behavior of (. The 
fact that they are regularized will be enough to show that they cannot induce a time dependence on 
C At this point we can also notice that has a tadpole in de-Sitter space which renormalizes the 



^"'^Noticc that wc can even assume that these quantities have been renormalized: contrary to the case where we 
integrate first in time up to plus infinity and then in momentum, here we are doing the integration in opposite 
order, and therefore all UV divergencies are true UV divergencies of the theory and have to be reabsorbable either 
by a counterterm or by a field redefinition. 
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background Universe and its contribution to is cancelled as explained in sec. 4.0.1. We can therefore 
concentrate only on the connected part of the correlation {{a'^)k{t){a^)k'{t'))- 
We can rewrite (((T^)fe(t)((T^)fe/(t')) as 



{{&'')t{f){&'')kif)) = / A/ <iVe'<"+'='*l((a=)(f, (')) (112) 

^ .«(.>J')/A_e---((*')(x-_.*,(.')(0.O). 

where x_ = |x — x'\/2. Since we are interested in small /c's if this integral converges even for A; = then 
we can set the oscillating exponential to one, and do the integral. In this way we will obtain a result 
which is exact up to corrections of order k/{a{t*)H) <C 1. We are therefore led to evaluate 

{i&\Mt)ia\,^oit')) = 6(^\k + k') J d3x_((a2)(f_,t)(c72)(0,t')) ■ (113) 

In order to do this, we take the correlation function {{a){x-,t){a'^){0, t')) which can be constructed from 
{a{x-,t)a{0,t')), and take the large distance limit. The resulting expression is 

{a{x-,t)a{0,t')) - H'^log{AiRAx^) , HAx > 1 . (114) 

where Ax^ = —(?? — ??')^ + ^~ A//j is the IR cutoff. Notice that at equal times this is the standard 
logarithmic dependence of massless scalar fields in de-Sitter space. If we define 

Iit,t') = J d^x.{ia^)ix.,t)icr')iO,t')) , (115) 

then we have that at large distances 

(d2(f_,t)(a2)(0,i))--7^ (116) 

which tells us that I{t,t) is convergent in the infrared. Since the expression is regularized (or even 
renormalized), then it converges even in the UV. This implies that I{t,t) is finite. Notice that in order 
for I{t, t) to be compatible with the symmetry 

a— >-Ao, x^x/X, k^Xk, (117) 



Ht,t)-^^. (118) 



it must be proportional to l/a{t)^. Dimensional analysis fixes the powers of H, and we conclude that 
I{t, t') goes as 

where cjjv is a number (expected to be or order one) that depends on the regularization and the renor- 
malization procedure, and that controls what is the strength of the correlations of 

Eq. (118) tells us that the correlation function of &^ is Poisson like distributed, i.e. its power spectrum 
is k independent. This is enough for us to be able to show that the contribution to ( from the cut-in- 
the- middle diagram is time independent and therefore scale invariant. 
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In fact by using that 

\l{t,t')\<V\I{t,t)I{t',t')\ (119) 
we can write the contribution of this term to eq. (110) as follows 



'^^'Hfc + ^')^e-3^*\ (120) 



where in the last passage we have neglected numerical factors. Eq. (90) implies 



t*=hog(-^] , (121) 



H \Heout 

so we obtain: ^ 

<5^'\k + k')^^cuvei,tN , (122) 

where we have reinserted the factor of N associated to the number of a fields. We see that the contribution 
from ((cj^)fc(t)((T^)fc/(t')) to the cut-in-the-middlc diagrams is time-independent and scale invariant. 

It is worth giving some physical interpretation of this result. It is quite easy to understand why the 
result is time-independent. In order for the CIM diagram to be able to induce a time dependence on (k, it 
is necessary for the self correlation of T(^cr)ixv to be coherent on a wavelength of order a{t) jk as the external 
time is taken to infinity. However, since the a fields are massless, the relevant fluctuations are only at 
most of order Hubble in size, and therefore Tf^a)^^ becomes quickly uncorrelated as the mode k spans 
many Hubble regions. This is why the main effect from the CIM diagram is peaked at the smallest time 
t* when the mode k spans the smallest number of independent Hubble patches. Notice in fact how the 
signal scales like e^^^ which is approximately the inverse of the number of independent Hubble patches in 
a box of radius a{t*)/k. Once the correlation function is time independent the symmetries of the problem 
forces it to be scale invariant. If we go back to real space, we have 

{U^,tf)ciM < ^^^uvNel^^j d'k p (123) 

which is invariant under rescaling of k. 

It is easy to see that the contributions from all the other terms in eq. (110) behave in a very similar 
way, and results in a time-independent and scale-invariant correlation function. We compute explicitly 
the contribution from the most interesting additional terms in App. B. We conclude that the contribution 
from the CIM diagrams are time-independent and scale-invariant and we proceed to study of the cut- 
in-the-side diagrams. 
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5.4 Late time contribution: cut-in- the-side diagrams 

We now proceed to the study of the contribution from the CIS diagrams. These diagrams require to 
compute the expectation value of the operator stress-energy tensor T(ct^^.^^ of the cr's due to the interaction 
with a C fluctuation at some earlier time. Let us represent this as 

(r(,)^,(x,t))^, . (124) 

Due to the divergencies involved in this calculation the expression of T^a)(;;iJ.u needs to be regularized. 
As in the former subsection we think of T(^(T)fj.u directly as the regularized, or even renormalized, stress 
tensor at that time. The one-loop CIS diagram is represented in Fig. 4 where the vertex labeled by 
t2 represents the perturbation to a cr mode due to the interaction with the primordial ( mode. There 
are three important things to notice that will make the proof possible even without having to explicitly 
perform the computation of the expectation value in (124). The first is that, because of translation 
invariance, the wavenumber k of the ( that perturbs the operator {Ti^cr)tii/)Ck ™ust be the same as the 
external ( for which we are computing the expectation value, and therefore it is very outside the horizon 
for the times of interest here: t > t*. The second important thing to notice is that because the cr's appear 
always with derivatives the a modes contributing to the finite part of {T(^^-^^y{x,t))(^^ have a momentum 
of order H or larger. The third important point is that at this order in perturbation theory the mode 
(" that perturbs T'(o--)^^(a;, t) is the free field that lives in the unperturbed metric. This implies that the 
initial Cfc approaches a constant in time value as k/{a{t)H) and that the induced metric approaches 
the unperturbed one, up to the simple rescaling of the scale factor 

a{t) e'^'=a(t) . (125) 

The above points can be summarized by saying that we are left to compute how the expectation value 
of {T(^„^^y{x, t)) gets altered in the presence of a metric that is just the unperturbed one with the simple 
rescaling of the scale factor, which is unobservable in the limit that it is spatially constant. We write 
'^{(j)ixv using the obervable quantities p, v in the standard way as: 

= ((Pa + P'^)^^'^^^ + P'^^i)Ck (126) 

where 



and 



{p.) = {pAt))o + {5pa)cu , (128) 

{Pa) = {Pa{t))Q + {Spa)<:k 

where the subscript ^ represent that these are expectation values taken with respect to the a fields in 
the background of ^, while the subscript o represents that they are computed in an unperturbed metric. 
The fact that a super-horizon wavelength becomes locally unobservable implies that for k/ {a{t)H) <C 1: 

(*p»(t))& ~ (ipAt))a ~ o (w(*)>o^&(!)) , (129) 
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There is one subtlety in this argument that we would like to put into evidence: when we compute 
{T(^a)tiu{x,t)), the insertion of the Cfc perturbation as represented in Fig. 4 can be at such an early time 
t2 <C t* that the mode ("^ is inside the horizon. In this case, the mode (" is observable, and in principle 
could imprint an effect on the T(^a)fiu- However, in this regime the modes whose frequency is higher than 
order at i ~ (which are the modes contributing to T(^cT)iJ,iy) much more energetic than the one of 
and therefore they just follow their adiabatic interacting vacuum up to t ~ A related subtlety 

in the above equation is also the fact that we are declaring that (5/9o-(i))Cfe) {^PfT{'^))Ck ' ^^'^ (<^^o-(0)Cfc 
depend on (" evaluated at the same t. This is of course not true, and in general the relationship will be 
non-local in time but given that at late times Ck{t) approaches a constant and that the expectation value 
is dominated by this regime we can approximate the dependence as being Ck{t)- 

It is easy to see that eq. (129) is enough for proving that the correlation function of ( becomes time- 
independent out of the horizon. In order to do this, we write again the continuity equation for the full 
stress tensor 

(VMr'^°)a = (130) 



which gives 



C+|3ff-2| + f C^C (131) 



+ ^ + ^ {{5p.)a+SHi{5p.),, + (Sp.),,)) 

+ 2^((Pa)o + (p.)o)(^H(i))a + ( 6i/2 + j 6N^,2 + Hxc.2 + HSk^,2 = 
where we have defined: 

N=1 + SN^ + 5N^2 + +5N(^^2 , (132) 
X = diN' = xc + + xc.2 , (133) 

with 6Ni^ being proportional to ^, 6N^2 to and 6N^^2 to (^a^, and similarly for x- Notice that in 
eq. (131), SN^2 and Xa^ do not appear. This is so because the we are taking the expectation value with 
respect to the a field, and the tadpole cancellation ensures that those terms cancel. 

As usual, we arc interested in cq. (131) when the mode k is well outside the horizon. In this regime, 
we can find solutions for dN^^2 and x^a^ at leading order in k/{a{t)H) <^ 1. By solving at this order the 



^^Noticc that this same kind of reasoning applies to contributions to the C correlation function at times after 
inflation, for example during the period of reheating or during a GUT phase transitions (if this exists). These 
are epochs during which we know very little of what is going on. Still, the time scale associated to the thermal 
fluctuations for the case of reheating, or to bubble collisions in a GUT phase transition, is expected to be much 
faster than the Cfc mode when this is inside the horizon. Obviously, these fluctuations stay in their adiabatic vacuum 
at those times. 



38 



(0, 0) and the (0, i) Einstein equations we find: 



-^^C-^ ^ -T^ Txi^Pak, (134) 



3i?2 + ij) 



,C1. 



Eq. (129) implies that 5N(^„2 is of order k/{a{t)H)C,k and XCa'^ is of higher order. Therefore, by looking 
back at equation (131), in the long wavelength limit, we can see that the sources for Q at late time decay, 
which is enough to see that the correlation function of C, will not depend on time and is scale invariant. 
In fact by taking the leading terms in k/{aH) ^ 1 and in slow-roll parameters in eq. (131), and using 
(129), we obtain: 

Ck{t) - fdt'^ (l - ^) T^H%{t')-^ ~ 7|^a(*) fdt'^ (135) 
Jt* H \ a{ty ) M|ji7 a{t')H M^^H Jt* a(t')H 

^ , . k \ _ Iff* ^ , , 

j-Ck{t) — —e ^* ~ 77^a(*)w , 



Here for simplicity we have dropped higher order corrections in the slow roll parameters. The result can 
be trivially extended to include those as well. Corrections suppressed by higher derivatives in eq. (129) 
and (131) are down by further powers of eout ^ 1- For the cut-in-the-side diagrams we therefore obtain: 



{C2,kmAt))cis ~ 5^^\k + k')l,-^Ne. 



A;3e2M4j^'-"* ' ^^^^^ 

where we have reinserted the number N of spectator fields. This is a time-independent and scale-invariant 
result. 

So far in this section we have neglected almost completely the non IPX diagrams that we studied 
in sec. 4.0.1. The tadpole cancellation procedure applies also in this case where we perform the mo- 
mentum integration first and the time integration second. Therefore, they are automatically zero once 
the background has been redefined according to (86). All the CIS diagrams give a time-independent 
scale-invariant result. 

Finally, we need to study the case where a C, generated by the first term of the right-hand side of (91), 
that we called Ci, correlates with a term generated by the second term on the right-hand side of (91), 
that we called ^2- By using the standard property of correlation functions, we have 

m,i{t)Ck'm\ ^ {{Ck,imk',i{t))Y" {{Ck^cwmY" ■ (m 

Since each of the terms inside the square root is time-independent and scale-invariant, we conclude that 
so is for the correlation between ("i and C2- 

This last result completes the analysis of all the contributions to the correlation of ( at late time and 
therefore we conclude that the result is scale-invariant and time-independent. The procedure to compute 
the loops that we used in sec. 4 is justified in yet another way. 
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5.5 Slow-roll corrections 



So far in the above calculation we have not included all the slow roll corrections, for example those ones 
that come from the wavefunction of the fields as the spacetime deviates from de-Sitter. It is rather simple 
to include the first corrections due to the slow roll parameters as follows. In the former subsection 5, 
we have seen that C^. at some late time is the result of three distinct contributions: the free evolution of 
(^k from some initial condition at time (which is in practice just a constant), the contribution which 
comes from the variance of T^^^^^j^ (which is exponentially peaked at ti,), and finally the contribution from 
the perturbation to the expectation value of T(^a)nv due to a primordial C, fluctuation (which is again 
exponentially peaked at t^). 

The last two of these contributions are exponentially sensitive only to the time t*, and so, in the case 
of some time dependence of i7, H,. . ., the result will be the same as the one we obtained, just replacing 
H, H,... with H{U), H{U), .... 

The same is true also for the first contribution if is chosen to be close to the time of horizon 
crossing so that the variation of the parameters of the spacetime between and the time of horizon 
crossing is neghgible. This is clearly possible given the slow variation of the parameters with respect to 
the exponential stretching of the modes. We therefore conclude that with the replacement 

H, H,... ^ H{Q, H{U),... , (138) 

we can incorporate the leading slow roll corrections. 

5.6 TT^ large self- interactions 

The above proof of the constancy of ( was done for the case of N spectator massless scalar fields. The 
same result holds for the case of the tt"^ interactions. Considering for example the cubic interactions, 
everything in the above proof proceeds in the same way, apart from trivial combinatoric factors, if, of 
the three tt's in the vertex tt^, we consider one to be the infiaton and the other two to be two spectator 
fields. In fact, in the above proof, there were only two crucial points. The first is that for the CIM 
diagrams,the correlation of ct^ and the similar operators was not scale-invariant, but rather Poisson like. 
There is no difference in this case (notice that because here we could ignore metric fluctuations the proof 
would be even simpler). The second important point of the above proof was that, in the CIS diagrams, 
{^Pa)c,k^ {^Po-)Ck^ ■ ■ ■ were all going to zero as k/{a{t)H) 0. Clearly this will be the same also in this 
case, as the free mode n'^ becomes a constant out of the horizon and vr is only derivatively coupled. It is 
again correct to assume that (Tr^)^^., , • • ■ are all going to zero as k/{a{t)H) 0. We therefore conclude 
that the correlation function of C, is constant also in this case, justifying also in this alternative way the 
calculation of sec. 3. Slow roll corrections can be included as before. 

6 Other interactions and higher loops 

Finally we comment on how our results are expected to generalize to higher loop calculations and for 
more generic interactions. 

The two kind of theories we have studied resulted in logarithms of the form log(H/ fi), where fi is the 
renormalization scale. Current bounds of the non-Gaussianities of the CMB that limit the size of the 
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TT^ self-interactions [14] make these corrections extremely small. The one loop effects change with the 
A;-mode only slowly due to the slight dependence of the Hubble scale on the horizon-crossing time of the 
A;-mode. These slow-roll suppressed effects are in general very small but it is possible to imagine that if 
this were not to be the case, then a procedure similar to the renormalization group could be developed 
and implemented. This lies beyond the scope of the current paper. 

One common characteristic of the two interactions we had studied was that the fields running in 
the loop appeared in the Lagrangian with at least one derivative. What happens when this is not the 
case, for example if the a fields have a (positive) mass term? This case has been studied and solved 
by Weinberg in [22]. In his first paper [2], he noticed that if the mass term was treated perturbatively 
as a two-line vertex, and one used in the loops the wavefunctions of a massless scalar field in de-Sitter 
space, then the integral in time in the loops would not converge. This interaction might then give rise to 
a time-depedencc of the C, correlation function. In a following paper [22], however, he solved the puzzle 
by noticing that summing up all the perturbative mass insertions into the propagator (i.e. using the 
wavefunction of a massive field in de-Sitter) resulted in a wavefunction that decayed exponentially at 
late times and the time integrals in the loop converged again. The apparent time-dependence of the C, 
correlation function disappeared. 

We would like to see how this effect appears in the arguments we presented in the previous section. 
For massive cr's there is a term proportional to m^cr^ in the stress energy tensor. Let us imagine treating 
the mass term perturbatively and repeating the steps of sec. 5. We expect to find a time dependence in 
the C correlation function. 

In the cut-in-the-side {CIS) diagrams of sec. 5.4, we would have to compute the perturbation to the 
expectation value of the a stress tensor T(^a)^iu due to a background free C, mode. It is still true also in 
this case that a constant (" is unobservable and that therefore the perturbations to and would be 
proportional to the first derivative of This is enough to make the source in the CIS diagrams shut 
down at late time, and therefore there is no induced time dependence from the CIS diagrams in the case 
of massive cr's either. 

The cut-in-the-middle {CIM) diagrams of sec. 5.3 are different. Those diagrams involved the cor- 
relation function of T(g.)^j^. In the massless case, the various terms of the stress tensor involved always 
derivatives acting on the a fields, and this was enough to ensure that the correlation functions of the 
various terms of the stress tensor became Poisson like for distances longer than the horizon. This implied 
that these fluctuations were not able to source a C mode once the wavelength was much outside of the 
horizon. In the massive case things are different. The part of the stress tensor that is proportional to 

does not involve derivatives of a which means that the correlation function of this part of the stress 
tensor is not Poisson but scale invariant. This part of T(^cr)ni/ has correlations for very long distances 
and therefore it is able to source a ^ mode regardless of how much outside of the horizon that mode is. 
Physically this is quite easy to imderstand. Let us consider a (" fluctuation of wavenumbcr k. At the time 
when this fluctuation crosses the horizon there will be a fluctuations of T^aji^u on the same wavelength 
generated by two a fluctuations one of which has a wavelength of order 1/k, while the other has a longer 
wavelength and crossed the horizon at an earlier time. At this point if we are using the massless a wave- 
functions the cr's freeze and remain constant. They remain up in their potential wherever they happened 
to land at horizon crossing. This results in a region of additional potential energy which now begins to 
expand as the Universe expands, without never decaying very much like a space-dependent cosmological 
constant where the space dependence is given by the wavenumber k of the T(^a}tJ.u fluctuation. Since the 
additional potential energy never decays it is able to coherently source the mode (k and therefore it 
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Figure 5: Mass insertions to obtain the massive propagator. 



induces a time-dcpcndcnt effect. Notice tliat ffuctuations of T(^a)i^iv the same scale produced at a later 
time by shorter scale ct's are instead suppressed, because they average out as in the massless case. It is 
just the fluctuations of the a modes of wavelength up to order k that matter. 

It is clear that as a result of the approximations used the above calculation does not capture what will 
actually happen. Once the mode is outside of the horizon the a field cannot remain up its potential and 
will move down and reach the minimum. In fact it will do so either exponentially fast in time undergoing 
oscillations (if m > 3H/2) or very slowly (if m < 3H/2). We just do not see this effect above because 
we have treated the mass term perturbatively. Since the out of the horizon solution for a massive field 
scales as [22]: 

a(i)^± , (139) 

where 



2 

the effect of the mass term is perturbatively of order Ht. If one is interested in times when this 
correction is important then the only way to properly take it into account is to re-sum all the mass 
insertions and use the wavefunctions of a massive field (sec fig. 5). If one were to do so, then, after some 
long time in the case of small m one would see that the source for the C mode shuts down and ^ becomes 
constant. Notice that for small masses this happens after a time 

Wether this time scale is long or not depends on the problem. For example, in standard inflation, this 
time scale might be longer than the time to reheating and in this case the time dependence of C will be 
relevant. Notice however that the overall effect goes down as the energy density, as mP, and therefore it 
becomes smaller and smaller in the limit m ^ 0. This time scale tc is clearly short if one wishes to study 
time-dependence of correlation functions in eternal inflation or even in de-Sitter (which would require 
changing gauge and defining some relevant observable [18]). In this case one can conclude that the ( 
correlation function becomes constant in time a time of order tc after horizon crossing. 

What happens if the a fields do not have a mass term, but for example a potential of the form 
or even higher dimensional potentials such as cr^/A^, where the cr = point is a minimum? All these 
terms would enter if we were to do the calculation at more than one loop. Let us just concentrate on 
Xa^, the other cases being just a trivial higher-loop generalization of this one. Clearly if one were to treat 
the A term perturbatively one could exactly reproduce the argument we just made for the massive case 
(extended this time to two loops) and find that the CIM diagram produce a time dependence in the ( 
correlation function. The solution is the same as in the massive case. In that case, we had to re-sum all 
the insertion of the mass term in the propagator (see fig. 5), here we will have to re-sum all the insertion 
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Figure 6: One-loop insertions to obtain the effectively massive propagator. 



of the quartic interaction in the propagator (sec fig. 6). This will generate effectively a mass term for 
the a fields and will make the field decay and the induced time dependence of the C correlation function 
disappear eventually. As in the former case it is still true that the time at which a (^^ mode becomes 
constant after horizon crossing might be very long. There is just one subtlety in the above argument. 
This is the fact that in principle a massive counter term could cancel the mass term induced by the 
quartic interaction (this would require a large amount of tuning, as we know from the Higgs particle in 
the Standard Model but still in principle it could happen) and therefore one would still have a massless 
field which would therefore induce a time dependence on Q by jumping up its potential and staying there. 
We find it very hard to imagine that this cancellation is actually possibly in a quasi de-Sitter space as we 
expect that finite terms in the loop to inherit some time dependence from the time dependence of H and 
the counterterm should not be able to reabsorb that. However even assuming that the mass can be tuned 
away there are other diagrams that are relevant and that make the field run down. In fact a particle in 
a quartic potential rolls down to the bottom even in the absence of any explicitly quadratic term. For 
example one could take the diagram in Fig. 7 where all cr's in the dashed lines which represent primordial 
fluctuations get correlated with other primordial ones from the right-handed part of the CIM diagram 
This diagram will let the a field roll down its potential and therefore make C,^ constant eventually though 
possibly a long time after it has crossed the horizon. 

Notice that similarly to the apparent time-dependence of ^ induced by these kinds of interactions, 
we expect there to be also IR divergencies associated with the momenta loop integrals. These have been 
found for example in [24, 27, 30, 31, 32] in the context of spectator scalar fields and of the standard 
slow rolling inflaton. If we consider for example the case of the massive spectator a fields we discussed 
at the beginning of the section, these IR divergencies are associated to the counting of all the a modes 
that crossed the horizon before the mode k did. If one uses the massless wavefunctions, these modes that 
crossed the horizon arbitrarily in the past do not decay and all contribute to the correlation function, 
inducing an IR divergency. As for the induced time-dependence of C, we expect that the resummation of 
the diagrams we discussed will let all the spectator fields decay, effectively multiplying the wavefunction of 
the modes outside of the horizon by some power of (A;r?), where k is the wavenumbcr of the mode. In this 
way, the IR divergency is expected to disappear as well. Reference [31] offers an example, valid in some 
particular theories, of the resummation techniques of the diagrams and shows that upon resummation of 
the wavefunction this particular kind of IR divergencies disappear. 

There are other interactions that Weinberg pointed out would make the time integrals in the loop not 

^^This diagram can be thought of as representing the difference between aa^ and cr((T^), where it is only the 
second term that can be potentially removed from the mass counterterm. 
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Figure 7: Left-handed part of one of the CIM diagram whose resummation will make a decay out of the 
horizon even in the absence of a mass term. Each of the cr's in the dashed lines has to be correlated with 
a corresponding one on the right-handed side of the diagram. 



converge and therefore apparently make C time-dependent . Studying these interactions goes beyond 
the scope of the current paper. We just mention here that an example of a term of this kind is (^a" with n 
being an integer. Such a term would for example arise when the a fields have a potential of the form cr". 
But because the cr's have a potential they would roll down to their minimum and would not induce a time 
dependence of (" even after including the ("cr" vertex. We believe this is missed in Weinberg's theorem 
[2, 22] as in its proof every interaction is analyzed on its own assuming that all fields are massless. In 
reality, it seems to us these kind of interactions are highly constrained by symmetries and come in groups. 
Therefore we think they should be studied together in specific setups. 

Though we believe the above arguments to be very sound, strictly speaking we did not prove them. 
It would be interesting to directly verify them and also to explicitly find and study cases, if they exist in 
local theories, where the kind of interactions we discussed in the last paragraph appear and do induce a 
time dependent effect of the C correlation function even after the proper treatment of the mass term and 
of the related diagrams we discussed in Fig.s 5, 6, and 7 has been taken into account. 



7 Conclusions 

In this paper we have studied one-loop corrections to the two-point function of the curvature perturbation 
C, in two kind of theories: one in which the inflation has large self- interactions, and one in which the 
inflaton interacts gravitationally with N massless spectator scalar fields. 

In both cases, we found that the one-loop corrections are time- independent and scale-invariant. Fur- 
ther, we have shown that there is a logarithm of the form 

(Cl)i-ioop~pX/31og(^^) , (142) 

where H is the Hubble scale during inflation, is the renormalization scale, and /3 can be thought of 

^''These are for example the ones that contain time derivatives of the fields, and that appear in the Lagrangian 
as multiplying in three dimensions a number of scale factors a equal to one, after counting as minus two powers of 
o every appearance of a time derivative [2]. 
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as being related to the /3 function of the correlation function. This result is physically sensible: loop 
corrections for a weakly coupled theory are small if we choose the renormalization scale close to the 
energy scale of the process we are interested in. For inflation, this energy is Hubble. 

Our results differ from the previous literature, which instead found that a logarithmic running of the 
form log(fe//i). After arguing that symmetries forbid this result we were able to identify the mistake 
in previous calculations. Dimensional regularization was applied incorrectly. We also performed our 
calculation using a different regularization scheme finding the same result. We were also able to explicitly 
renormalize the simplest of the two theories we studied. 

For the first time we studied the time dependence of the correlation function finding that for 
the interactions we considered C becomes constant outside of the horizon even at one-loop level. This 
is the case once the background spacetime has been properly redefined in order to take into account 
tadople diagrams. These findings offer an alternative proof that the divergencies have to be reabsorbed 
by local counterterms and are not instead infrared divergencies that appear as UV divergencies because 
of the particular way the calculation is performed. We achieved this not only by explicitly renormalizing 
the theory but also by performing the loop calculation at finite external time. We further developed 
an alternative method for computing loop corrections based on solving perturbatively the Heisemberg 
equations for the operators, which allowed us to work at every time directly with the renormalized physical 
quantities. In this approach the physical reason why ( becomes constant outside of the horizon appears 
clearly. It is due to the combined effect of the fact that a constant C mode is unobservable and that the 
stress energy tensor which sources ^ becomes rapidly uncorrelated for scales longer than the horizon. 

We are then able to use the intuition we developed to infer how our results generalize to different 
interactions and to higher loops. In particular, contrary to claims in the literature, we argue that no long 
term time dependence is generated by any form of spectator field, even if they have a potential term. 

Wc find our result to be important not only because generally they increase our understanding of 
infiation and of de-Sitter space, but also specifically for two problems. The first is related to the fact that 
during the evolution of the Universe from the end of inflation up to now there are various epochs, such as 
reheating, where we do not know the details of the physical processes at play and thus we could imagine 
that there were very large fluctuations on the horizon scale. The theory of the spectators massless scalar 
flelds can be thought of as a weakly coupled version of this, and it would have been worrisome if we 
had found a time-dependence on ( on large scales. In fact, while the effect from the a flelds would have 
been in any event small, the effect from fluctuations during reheating could have altered the value of ( 
completely and the predictivity of inflation would have been lost. Our flndings prove that this is not the 
case. 

A second reason why we consider our flndings relevant is connected to slow roll eternal inflation. If 
we trivially extend our finding to the inflaton two-point function at coincidence, the fact that the one 
loop corrections arc tinic-indcpcndcnt and scale-invariant means that the smoothed two point function 
of the infiaton field at coincidence grows linearly with time. This means that the results of [18, 19] still 
apply after including the kind of loop corrections we computed, and that therefore, at least for what 
these interactions arc concerned, standard slow roll inflation exists and there is a sharp phase transition 
as a function of the slope of the potential where eternal infiation begins. 

We think our results motivate other interesting calculations. For example, it would be interesting to 
verify with an explicit calculation our arguments of sec. 6 about the fact that the resummation of higher 
loops makes the time dependence induced by spectator fields with a potential at a fixed perturbative 
order disappear, following and generalizing the calculation of [31], and understand the issue also in the 
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case in which the inflaton and the graviton are let run in the loops. This would be particularly important 
for assessing the existence and the characterization of slow roll eternal inflation. Further, it would be 
interesting to study loop corrections and how the renormalization procedure works for the most generic 
theories of inflation. For single field inflation, this would amount to study the most general form of 
the Lagrangian presented in [3], of which we studied only one particular limit in sec. 3. There are 
also questions at a more fundamental level that our calculation might help answer. For example, in the 
dS/CFT correspondence [33, 34], it is speculated that the time-dependence of correlation functions, which 
occurs at pcrturbativc level in theories where spectators scalar fields have a potential term, should be 
mapped into renormalization group running for some operators in the CFT. It would be very interesting 
to see how this happens in a concrete way. Further, it is possible that the structure of the loop corrections 
that we found might be related, in string theory, to the density of single-string states, similarly to what 
happens in negatively curved spaces [35]. This is another direction worth exploring. 
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A Loop integrals with momentum frequency cutoff 

Here we give some details of the calculation of the loop integrals that lead to eq. (34). The two momenta 
integrals can be done exploiting the identity: 

/ ^ 2-71 /"^"(^l) fki+k 

(fkid^k2d'^^\ki + +k2 + k)f{h,k2,k) = — dkiki dk2k2f{ki,k2,k). (143) 

k Jo J\ki-k\ 

The result of the two momenta integrations therefore splits into the sum of two terms: one in which we 
take the contributions of small ki's such that ki — k < 0, and one where we take the contribution from 
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the large fci's such that ki — k > 0. The contribution from the large /ci's gives: 



X |i? 



,2 



drj 



4 4 
^1 



e^'^^'i (8 - llA;(r/i - r/2) + Ik'irH - r?i)) 



+ 



oAf !Z2 _i 

+e ^Ui 



^e^^'^i {h^tiI {hSI + 2Ai7r/i(r/i - 172) + 2A2(r7i - rya)^) (s - 5(?7i - 772) + k" {m - 772)')) + 

+4A3i?77i(77i - 772)3 (2 - A;(77i - r/2)) + 2A4(r?i - 772)^) + 
+e2^''2 (^H%1 {hWi + 2Aii-77i(77i - 772) + 2A^{rii - mf) (s + 5(771 - 772) + (r?i - 772)')) + 
+4A3i?r7i(r7i - 772)^ (2 + k{r]i - m)) + 2A\m - m)^)] } ■ 

(144) 

It is quite straightforward to realize that integrations of the terms in the third, fourth, fifth and sixth 
lines of (144) can not give rise to any logarithmic diverg ency. Because of the Exp (2^ - l)) the 
integrand becomes exponentially small as soon as I771 — 772 1 becomes larger than r]2H/A, which is the 
minimum distance in conformal time allowed by our cutoff. This means that for these term the integral 
is peaked at the cutoff region, and it can not give rise to any logarithmic dependence. Integrations of 
the terms in the second and third line of (144) instead is supported on all scales, and a straightforward 
integration leads to our result in (34). 

We do not reproduce here the contribution from the small fci's. It is quite straightforward and it 
gives rise only to finite terms. 

B Additional correlation functions in the CIM diagrams 

Here we complete the calculation of the correlation functions in (110) for the CIM diagrams. It is easy to 

see that the contribution to the ( correlation from (((j^)fc(*) {{diO')'^) (t')) and from (((Sjcr)^)^ (t) ((^jcr)^)^, {t')) 

proceeds in very similar terms to the term studied in (122), and results in a contribution parametrically 
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equal. To analize the contribution from ((gj^* {adia))^ (t) (^c^* (<^i9i(7))^, {t')) we have to study 

< 5^^\k + k') j d'x^ {{adia) (adia) iO,t)) 

X [-l^J ^(^^^-^^ {x-,t'){adma){0,t')) 

^S^^\k + k') j (fx- {&ix-,t)&iO,t)){dia{x-,t)dia{0,t)) 

p- / d?x- {a{x-,t')a{0,t')){dja{x-,t')dma{0,t')) 



where we have used (114) and in particular that {a{x-,t)dia{0,t)) = 0. Plugging back in the expression 
for {CO CI M, and proceeding as in the former example, we obtain that this term contributes in the 
following way 

{C2,k{t)C2,k'{t))ciM < S^'Hk + k'f-"^^ fdh f'dh ^ 



We find that this contribution is also time independent and scale invariant. Proceeding analogously with 
the correlation from {{o'^)^. (t) (gi^* {(Tdia))^ {t')) we find that 



which leads to 



fe.WfeH*)>c,M<*''H* + *022^^;*/;*;;(jp7^ (147) 



1 



which is again time-independent and scale-invariant. It is straightforward to see that the contribution 

from (^(^icj)^^^ {t) {-^9^ (adia))^, {t')) gives a parametrically equal contribution. This last one exhausted 

all contributions to the late time correlation of ^ from the CIM diagrams, and we therefore conclude 
that this contribution is time-independent and scale-invariant. 
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